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CV' ABSTRACT 

O ' 

o : 

We present the results of an investigation of the effects of Far Ultraviolet 
Q_i! (FUV) radiation (6.0eV < hv < 13.6eV) from hot early type OB stars on clumps 

in star-forming molecular clouds. Clumps in FUV-illuminated regions (or pho- 
todissociation regions or PDRs) undergo external heating and photodissociation 
as they are exposed to the FUV field, resulting in a loss of cold, molecular clump 
mass as it is converted to warm atomic gas. The heating, if rapid, creates strong 
photoevaporative mass flows off the clump surfaces, and drives shocks into the 
clumps, compressing them to high densities. The clumps lose mass on relatively 
short timescales. The evolution of an individual clump is found to be sensitive 
^ | to three dimensionless parameters: r] c0 , the ratio of the initial column density of 

the clump to the column iVo ~ 10 21 cm -2 of a warm FUV- heated surface region; 
v, the ratio of the sound speed in the heated surface to that in the cold clump 
material; and tFuv/t c , the ratio of the "turn-on time" tpuv of the heating flux 
on a clump to its initial sound crossing-time t c . The evolution also depends on 
whether a confining interclump medium exists, or whether the interclump region 
has negligible pressure, as is the case for turbulence-generated clumps. In this 
paper, we use spherical 1-D numerical hydrodynamic models as well as approx- 
imate analytical models to study the dependence of clump photoevaporation on 
the physical parameters of the clump, and to derive the dynamical evolution, 
mass loss rates and photoevaporative timescales of a clump for a variety of as- 
trophysical situations. Turbulent clumps evolve so that their column densities 
are equal to a critical value determined by the local FUV field, and typically 
have short photoevaporation timescales, ~ 10 4-5 years for a 1 M Q clump in a 
typical star- forming region {t] c q — 10, v — 10). Clumps with insufficient magnetic 
pressure support, and in strong FUV fields may be driven to collapse by the com- 
pressional effect of converging shock waves. We also estimate the rocket effect 
on photoevaporating clumps and find that it is significant only for the smallest 
clumps, with sizes much less than the extent of the PDR itself. Clumps that are 
confined by an interclump medium may either get completely photoevaporated, 
or may preserve a shielded core with a warm, dissociated, protective shell that 
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absorbs the incident FUV flux. We compare our results with observations of some 
well-studied PDRs: the Orion Bar, M17SW, NGC 2023 and the Rosette Nebula. 
The data are consistent with both interpretations of clump origin, turbulence and 
pressure confinement, with a slight indication for favouring the turbulent model 
for clumps over pressure-confined clumps. 

Subject headings: ISMxlouds - ISM:dynamics - stars:early-type - stars:formation 
- ultraviolet: stars 



1. Introduction 

Young massive OB stars significantly influence the structure, dynamics, chemistry and 
thermal balance of their associated molecular clouds through the impact of their ultraviolet 
photons. Their extreme ultraviolet photons (EUV, his > 13.6 eV) ionize the gas immediately 
surrounding them, forming H II regions. Far ultraviolet photons (FUV, 6 eV < his < 13.6 
eV) dissociate molecular gas beyond the H II region, creating a photodissociation region or 
PDR (Hollenbach & Tielens 1999). PDRs are ubiquitous: FUV photons dominate the gas 
heating in PDRs and affect the chemistry and physics of gas over a large fraction of the 
volume and mass of molecular clouds. The study of interactions between FUV radiation and 
molecular cloud gas is therefore important in understanding molecular cloud evolution and 
the feedback between massive star formation and subsequent star formation in molecular 
clouds. 

Observations of PDRs compared with theoretical models probe the physical conditions 
in star-forming regions through dust and gas emission mainly in the infrared (IR) wavelength 
region of the spectrum. PDRs are the source of most of the IR emission in the Galaxy. Dust 
grains and large carbonaceous molecules such as polycyclic aromatic hydrocarbons (PAHs) 
absorb radiation from the stars, and re-radiate this energy flux in the infrared, producing 
a continuum with solid state and PAH spectral signatures. FUV photons incident on dust 
grains and PAHs also cause photoelectric ejection of electrons which then collisionally heat 
the gas, a process called the grain photoelectric heating mechanism (e.g., Watson 1972, Bakes 
& Tielens 1994, Weingartner & Draine 2001). The gas cools via emission in many IR lines. 
Emission from PDRs has been quite successfully modelled (see Hollenbach & Tielens 1999 
and references therein). The models are usually parameterized by Go, the ratio of the FUV 
flux to the Habing (1968) FUV band flux of 1.6 x 10~ 3 erg cm~ 2 s~ 1 characteristic of the 
local interstellar radiation field, and n, the density of the gas. For advecting (non-stationary) 
models, the flow velocity of the material through the PDR is also a parameter. The FUV 
fields in PDRs near OB star- forming regions such as the Orion Bar or Ml 7 SW are typically 
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very high, (Go ~ 10 4 ~ 5 ); the PDR gas in these regions is inferred to be very dense with 
average densities (n) ~ 10 4 ~ 5 cm -3 (Tielens & Hollenbach 1985b). On larger scales and in 
more evolved regions such as the Rosette Nebula, the FUV fields and the average densities 
may be significantly lower by factors of 10 2 — 10 3 . 

Infrared, sub-millimeter and millimeter wavelength observations of PDRs indicate that 
gas in PDRs is inhomogeneous (van der Werf et al. 1993, Luhman et al. 1998), which is 
not surprising since molecular clouds themselves are observed to be highly inhomogeneous 
on a wide range of scales from tens of parsecs down to tenths of a parsec (e.g., Genzel 1991). 
There is both indirect and direct evidence for dumpiness in PDRs. Spatially extended fine 
structure line emission of neutral and singly-ionized carbon near H II regions shows that the 
FUV penetrates deeper into the cloud than predicted by homogeneous models, suggesting 
that the gas is clumpy in nature (Meixner et al. 1992, Stutzki et al. 1988, Steiman-Cameron 
et al. 1997). Further, strengths of high excitation lines of molecular species like CO J=14 — >13 
are observed to be greater than that predicted to arise from gas with density inferred from 
other PDR species. It is thus inferred that in regions such as M17 SW and the Orion Bar, 
the PDR gas consists of a "low" density (typically ~ 10 3 ~ 5 cm~ 3 ) component responsible for 
the extended CI and CII emission and high-density clumps (n ~ 10 6 ~ 7 cm~ 3 ) which give 
rise to the high excitation CO lines from their surfaces or, more precisely, at a depth of 
A v ~1 from their surfaces (Burton et al. 1990, Meixner et al. 1992, Steiman-Cameron et al. 
1997, Hogerheijde et al. 1995). Clumps are also observed directly by mapping PDRs at high 
spatial resolution in the 2 fxm lines of H 2 (e.g., Orion Bar, Luhman et al. 1998; Eagle Nebula, 
Allen et al. 1999). The small regions of enhanced emission are interpreted as indicating the 
presence of high density clumps (n ~ 10 6-7 cm~ 3 ) with sizes < 0.02 pc, occupying a very 
small volume filling factor ~ 1 — 2%. However, this interpretation is not very conclusive, 
and in this paper we present evidence to show that small clumps with such high densities 
may not exist in PDRs. Small differences in temperature within the observed regions can 
also explain the observed emission, as H 2 emission is highly sensitive to gas temperature 
(Marconi et al. 1998). 

Dense PDRs like the Orion Bar are spatially thin structures and if these PDRs are 
indeed clumpy, the clumps must have sizes smaller than the thickness of the warm PDR, 
~ 0.03 pc for the Orion Bar. The presence of small clumps affects the infrared spectrum from 
PDRs by introducing a range of densities in the FUV-illuminated region and by introducing 
enhanced advection in clumps that are photoevaporating. One of the goals of this paper is 
to investigate the conditions under which small, high-density clumps exist in PDRs. 

Small clumps might not exist in significant numbers if the intense FUV field photoe- 
vaporates them quickly. The incident FUV field heats up and pressurizes the surface layers 
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of the clumps, causing the layers to expand. This mass loss may cause a complete photoe- 
vaporation of the dense clumps, and thus destroy them on short timescales. The mass loss 
timescales due to photoevaporation of an FUV-exposed clump can be estimated by assuming 
that the heated gas at the surface flows outwards at its sound speed, cpdr- The mass density 
at the base of this heated flow, pb, is lower than the mass density, p c , in the cold clump, as 
we will discuss below. (See Table 1 for a list of frequently used symbols in this paper.) The 
mass loss rate is given by dm c /dt ~ ^Trr^pbCpoR and since m c ~ p c r c; the photoevaporation 
timescale is tpp ~ m c / \dm c j dt) ~ r c p c / (c PDR p b ). For a typical clump, and an FUV-heated 
surface at 1000K, this is approximately 3 x 10 3 (r c /0.01pc)(p c /p b ) years. Asymmetric mass 
loss from the surfaces of very small clumps can rocket them (e.g., Oort & Spitzer 1955) 
out of the PDR back into the cloud. Clumps that survive the rocket effect and remain in 
the PDR are likely to lose significant fractions of their mass on short timescales of order 
< 10 4-5 years, depending on the ratio p c j pi>. The clump evaporation timescale can be longer 
or shorter than its residence time in the FUV-illuminated surface region of a GMC. From a 
frame of reference moving with the advancing ionization front, interclump material and the 
clumps in the GMC are advected into the PDR and ultimately flow into the H II region. 
Large clumps survive and they can penetrate the H II region and affect its evolution and 
structure, becoming Evaporating Gaseous Globules or "EGGs" (Bertoldi & McKee 1990). 

In the present paper, we focus on the heating of an individual clump in a PDR by 
FUV photons from a nearby O or B star. We study a range of FUV fluxes incident on a 
range of clump sizes and determine the photoevaporative lifetimes of the clumps and the 
evolution of their structure. Our investigation is analogous to the study by Bertoldi (1989) 
and Bertoldi and McKee (1990) on the effects of EUV heating of clumps in H II regions, and 
a generalization of the studies by Johnstone et al. (1998) and Storzer & Hollenbach (1999) 
of the evaporation of small clumps and protoplanetary disks ("proplyds", O'Dell et al. 1993) 
by EUV and FUV photons. 

In a subsequent paper (Gorti & Hollenbach, in preparation) we will discuss the cumula- 
tive effect of O and B stars on Giant Molecular Clouds (GMCs), and the relative importance 
of EUV and FUV photons in dissociating and destroying GMCs. Ultimately, we are in- 
terested in the role of FUV radiation in regulating star-formation in GMCs. Destruction 
through photodissociation and photoevaporation of potentially star- forming clumps/cores in 
GMCs reduces the star-forming efficiency of GMCs, and the FUV field may thereby play an 
important role in regulating or limiting star formation. Photoevaporation of the outer enve- 
lope of a collapsing isothermal sphere by FUV radiation could rapidly deplete the material in 
the envelope, and may limit the final mass of the star formed. On the other hand, compres- 
sion of non-collapsing clumps by shock waves driven by the warm surface gas could possibly 
drive the inner cores to instability and gravitational collapse, triggering star formation. 
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In §2 of this paper, we discuss the present understanding of GMC structure and existing 
models. In §3, we describe the evolution of an H II region and its surrounding PDR as they 
propagate into the GMC. Timescales relevant for the FUV heating of a clump in various 
contexts are briefly discussed in §4. In §5 we present a simple analytical model which provides 
an understanding of the underlying physics of FUV-induced clump evolution. We consider 
various possible evolutionary scenarios and their dependences on clump parameters and the 
strength of the FUV field. For the more general case, a 1-D numerical hydrodynamics code 
has been developed and these results are described (§6). In §7, we estimate the significance 
of the rocket effect due to asymmetric mass loss and its effect on clump lifetimes in the PDR, 
and in §8 we discuss the implications of our results for observations of PDRs in star-forming 
regions. Application to observed PDRs is made in §9. We conclude with a summary of the 
present investigation (§10) 
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Table 1: List of symbols frequently used in this paper 



Symbol 


Meaning 


c c 


Thermal sound speed in cold clump gas 


CPDR 


Sound speed of heated, warm PDR gas 


G 


Strength of the FUV field in Habing units, 1.6 x 10~ 3 erg s" 1 cm -2 


N 


Fiducial column density equal to 2 x 10 21 cm" 2 


n c o 


Initial gas number density of clump 


m c o 


Initial mass of clump 


r c o 


Initial radius of clump 


tc 


Sound crossing time in cold clump gas, r c0 /c c 


t e 


Expansion timescale for clumps with r] c0 > i] crit 


tFUV 


Timescale on which the cold clump gas gets heated to the maximum temperature 




in a given FUV field 


tpE 


Photoevaporation timescale 


t a 


Shock-compression timescale for clumps with i] c o < r\ cr n 


VlF 


Velocity of the Ionization Front moving into the GMC 


XpDR 


Thickness of the PDR in the interclump medium of the GMC, bounded by the IF and DF 


a 


Initial ratio of turbulent to thermal pressures in the cold clump gas 


P 


Initial ratio of magnetic to thermal pressures in the cold clump gas 


7 


Power law exponent of the variation of density with magnetic pressure 


So 


Initial thickness of PDR shell on the surface of a clump (= N /n c o) 


Vcr 


Critical column density for complete photoevaporation of pressure-confined clumps 


VcO 


Normalized initial column density to centre of clump, n c0 r c0 /iVo 


T)crit 


Critical column density in a given FUV field towards which unconfined clumps evolve 


Vf 


Final column density of pressure-confined clumps in an FUV field 


X 


Photoevaporation parameter, (1 + a + /3)(r] c o — l)/(2(2z/ 2 + a)) ; A = 1 for rj = f] crit 


V 


Ratio of sound speeds in warm and cold gas, cpdr/c c 
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2. Clumps in Giant Molecular Clouds 

Giant Molecular Clouds are observed to have a hierarchical substructure that is highly 
inhomogeneous, with density enhancements such as cores, clumps and filaments on all ob- 
servable scales, from ~ 100 pc down to ~ 0.1 pc. We use the term " massive cores" to denote 
large scale density enhancements (> 1 pc in size) and "clumps" to denote smaller scale struc- 
tures (~ few tenths of pc or less in size) 1 . GMCs and their more massive constituent cores 
are strongly self-gravitating, and are supported against collapse by turbulence and magnetic 
fields. The internal velocity dispersions of the substructures in molecular clouds are found to 
scale with their sizes, with decreasing non-thermal support on the smallest scales (Goodman 
et al. 1998). The column density through a GMC or a typical GMC substructure is found 
to be approximately a constant iV ~ 10 22 cm" 2 , although significant variations (N ~ 1 fj 21 " 23 
cm -2 ) about this typical column may occur. (For a recent review on GMCs, see McKee 
1999). 

The smaller substructures or clumps are not all gravitationally bound and the exact 
physical nature of these clumps remains uncertain (Williams et al. 2000). They have been 
interpreted as temporary fluctuations in density caused by supersonic turbulence within 
the cloud (e.g., Falgarone & Phillips 1990, Scalo 1990). This view is supported by their 
apparent fractal nature, which seems to show the same structure on smaller scales down to 
a few tenths of a parsec (Elmegreen 1999), and their observed supersonic linewidths (Blitz 
& Stark 1986). In this picture, the gravitationally unbound clumps are transient objects 
formed by converging supersonic turbulent flows. The non-gravitating clumps are therefore 
being constantly formed and destroyed, and have short lifetimes of the order of their sound 
crossing time, t c ~ t c0 /c c , where c c is the clump sound speed, and r c0 is the initial clump 
radius. The effect of the sudden turn-on of an FUV field on such clumps is to drive a 
shock into the clump, producing a smaller, denser clump with a warm evaporative outflow. 
Surprisingly, the evaporative timescale is now somewhat longer than t c , as we will show. 
Assuming that the FUV flux does not strongly affect the turbulence so that the formation 
timescale remains the same, the presence of an FUV field will not greatly change the non- 
gravitationally bound clump abundance, but it will compress them and alter their thermal 
structure and dynamical evolution. 

Another interpretation is that small bound clumps as well as small unbound clumps 



1 Note that our terminology of cloud substructure is the reverse of what is sometimes used in the literature. 
We call the large scale structures " massive cores" and the small scale structures "clumps" to be consistent 
with observational work on dense regions in PDRs, wherein the density enhancements are usually called 
"clumps" . Our notation is also consistent with the notion of "hot cores" , where massive stars form. 
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represent stable physical entities confined by interclump pressure. Observations through 
molecular line spatial maps and velocity channel maps or position- velocity diagrams indicate 
the presence of dense coherent structures. Indirect evidence suggests the presence of a low 
density (~ 1 — 2 orders of magnitude less than the clumps) "interclump medium" (ICM) 
(see Williams et al. 1995 and references therein). The ICM exhibits high velocity motions 
(Af pa 10 km s^ 1 ) and therefore may have significant pressure, which could be of thermal, 
magnetic or turbulent origin. Clumps in this picture are confined by the ICM, and some of 
them may be self-gravitating (Maloney 1988, Bertoldi & McKee 1992). Pressure-confined 
clumps live much longer and presumably form more slowly than the gravitationally unbound 
clumps of the turbulent model. Therefore, FUV heating of these clumps and their destruction 
through photoevaporation will have a much more pronounced effect on lowering the clump 
abundance in this interpretation. 



3. The Overall Evolution Of A Massive Star Forming Region 

Massive stars form in the central regions of clumpy, massive, cores (M > 1OOOM ), 
of radius ~ 1 pc inside GMCs (see review by Evans 1999). These massive cores contain 
clumps whose densities may reach n c ~ 10 5 ~ 7 cm~ 3 , or even higher (Plume et al 1997), and 
the clumps are surrounded by interclump gas with lower, but uncertain densities, possibly 
nicM ~ 10 3-5 cm -3 . When a massive star forms, the EUV photons from the star immedi- 
ately ionize the nearby interclump gas, and form an H II region. The clumps in the H II 
region are also exposed to the ionizing radiation and evaporate (Bertoldi 1989, Bertoldi and 
McKee 1990). The FUV photons penetrate further through the interclump gas beyond the 
interclump ionization front (IF) and dissociate the molecular gas, forming a dissociation 
front (DF) and a PDR. 

In the earliest stages of evolution when EUV photons emerge from the forming O/B 
star, the IF moves as an R-type front. However, the IF quickly stalls and becomes D-type, 
at which point the high pressure H II region drives a ~ 10 km s" 1 shock into the neutral 
interclump gas beyond the IF. It is instructive to obtain a sense of sizescale at the R/D type 
transition. Assuming the massive star emits & Lyman continuum photons per second, the 
Stromgren radius at this juncture is r s pa 0.1<f)\g 3 n e ^ 3 pc, where 049 = 0j/lO 49 s _1 and the 
electron number density n e4 = n e /10 4 cm -3 . If the central star has an FUV luminosity of 
Lfuv — 10 5 L Q (typical of stars with 49 ~ 1), the FUV flux incident on the surrounding PDR 
is characterized by Go ~ 10 5 . The thickness of the PDR is approximately 10 22 cmT 2 /tiicm 
or 0.3n7c M4 pc. The DF initially advances ahead of the shock, but as photodissociation 
approaches equilibrium with H 2 formation, the DF slows down. Eventually, the shock moves 
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ahead of the DF, and the entire PDR is contained in the post-shock shell. The shock now 
impacts cold molecular interclump gas. Hill & Hollenbach (1978) showed that at this stage 
the shock velocity has slowed to < 3km s _1 , and that, if Uicm ^ 10 4 cm -3 , the shock 
completely stalls in t < 10 5 years. Therefore, there is an initial stage (t < 10 5 years, for 
dense ICM) in the evolution of the clumpy PDRs surrounding H II regions where the clumps 
are both shocked and exposed to the photoevaporating effects of the FUV radiation. This 
can be followed by a stage where the shock has stalled and unshocked clumps move into the 
FUV-illuminated zone. 

Two configurations are possible after the shock stalls. If the H II region is still completely 
embedded, it becomes stationary with its pressure matched by the molecular cloud pressure 
and the ionizing photons completely absorbed by recombining electron/proton pairs in the 
H II region. Much more likely is the situation where the H II region breaks through the 
surface of the cloud and becomes a "blister" H II region. In this case, the ionized gas can 
expand away from the cloud, and into the interstellar medium. The IF then slowly eats its 
way into the PDR and the DF similarly advances into the molecular cloud. The velocity of 
the IF advancing into the PDR is of order vjf ~ 1 km s -1 , just enough so that the flux of 
particles through the IF can balance the heavily attenuated flux of EUV photons reaching 
the IF (see Whitworth 1979, Bertoldi & Draine 1996, Storzer & Hollenbach 1998). In such a 
case, the clumps entering the PDR may not have been previously shocked and compressed. 
Further, in the turbulent GMC model, clumps form and dissipate continually and those 
formed in the PDR after the passage of the shock are not affected by it. In both the PDR 
around the Trapezium stars in Orion and in the M17 SW PDR (Meixner et al. 1992), there 
is no evidence for velocity shifts between PDR gas and the ambient cloud gas, or, in other 
words, no evidence for a shock at the present time. In this paper, we focus only on the 
photoevaporation of clumps, whether or not they have been shocked, and do not follow the 
potential shredding or flattening of a clump by the passage of a shock wave (see, e.g. Klein 
et al. 1994). 

The thickness or column depth of a PDR has been defined by Tielens & Hollenbach 
(1985) to include all the gas and dust where FUV photons play a significant role in the gas 
heating or chemistry. By this definition, in regions like Orion or Ml 7, the PDR extends 
to columns of ~ 10 22 cm" 2 and includes gas which is cool and molecular (H 2 , CO), but 
where FUV still photodissociates 2 and H 2 0. We are interested in the surface column N 
of the PDR which is heated to the highest temperatures (T max ~ 100 — 3000K, depending 
on G and n). For G /n > 4 x 10~ 2 cm 3 , N « 1 - 3 x 10 21 cm" 2 (equivalent to a dust 
FUV optical depth of order unity) and is nearly equal to the column where atomic hydrogen 
converts to H 2 in stationary PDRs (Tielens & Hollenbach 1985). We shall call this region 
the "warm PDR" hereafter to distinguish it from the entire PDR. The thickness of the warm 
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PDR, No/n/cM, is of order 0.03 — 0.1n I ^, M4 pc. The crossing timescale for the IF through 
the warm PDR is tpDi? ~ 0.3 — 1 x W^n-jcMA^fb yfi^xs, where Vfs — vif/{\ km s 1 ). Gas 
and dust (and clumps) from the opaque molecular cloud interior are therefore advected on 
these timescales from the shielded cloud, through the warm PDR, and across the IF into the 
H II region. Figure 1 shows a schematic diagram of FUV radiation from nearby O /B stars 
impinging on massive cores (large sizescales, e.g., Rosette) and clumps (small scales, e.g., 
Orion). 
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Fig. 1. — A schematic diagram of PDRs on large scales (e.g. Rosette) and small scales (e.g. 
Orion Bar). The figure shows blister H II regions formed by the massive stars, and the PDR 
bounded by the ionization and dissociation fronts. Clumps in the PDR are exposed to FUV 
radiation and lose mass by photodissociation and/or photoevaporation. 
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4. Impulsive heating versus slow heating 

FUV-illuminated clumps are heated at their surfaces on timescales tpuy that could 
be either fast or slow compared to the internal sound crossing timescale within the clump, 
t c = r co/ c c- Impulsive heating implies that the FUV heating of the clump surface layer, 
characterized by column N , occurs with tpuv ^ t c . m this case, the outer region is heated 
more quickly than the entire clump can respond, and attains a much higher pressure than the 
cold central region. A shock propagates inward and a heated photoevaporative flow expands 
off the surface. 

Several situations in clumpy molecular clouds may lead to tpuv ^ t c . For a clump with 
r c0 ~ 10 17 cm and c c ~ 0.3 km s" 1 (T ~ 10K), t c ~ 10 5 years. The OB star may "turn-on" 
on a timescale < 10 5 years. Even after this initial flash of FUV flux, the clumps may shadow 
each other so that tp\jy is given by the time for a clump to move out of shadow, which is 
of order r c0 /v c , where v c is the clump velocity. In a turbulent, supersonic medium v c may 
be larger than c c so that tpuy < t c . Similarly, if clumps are formed in turbulent gas by 
supersonic converging flows, the FUV turn-on timescale is roughly the time for the clump to 
form, which again is tpuv ~ r c0 /v t , where v t is the supersonic turbulent speed. Therefore, 
impulsive heating is likely the most appropriate approximation for the turbulent model of 
clumps in molecular gas. 

Let us consider the other situation of long-lived clumps confined by the ICM. Here, slow 
heating is generally the best approximation, but in certain conditions, impulsive heating 
may be more appropriate. A common situation after the turn-on of the OB star may be the 
advancement of the ionization front and PDR front into the opaque molecular cloud, which 
occurs at speeds v IF ~ 1 km s -1 . In this situation t FUV is the time for the DF and the IF 
to move so that the clump is one FUV optical depth closer to the IF, or tpuv ~ tpon = 
Xpdr/vif, where Xpp>R = Nq/ujcm is the thickness of the heated PDR zone in the ICM 
of density tijcm- The criterion for impulsive heating tpjjy < t c , is more easily met by larger 
clumps that have longer sound crossing timescales, as tpuv is the same for all clumps. This 
gives us a limiting size for clumps likely to be impulsively heated, 

r c0 > X PDR ^ « (0.3 - 1)X PDR , (1) 

VIF 

which is of the order of the thickness of the interclump PDR region itself. Clumps with radii 
smaller than this are slowly heated as they emerge from the opaque molecular cloud into the 
FUV- irradiated regions. Since we only consider clumps within PDRs that have r c0 < X P p, R , 
pressure- confined clumps are generally slowly heated in our analysis. 

Clumps that undergo slow heating, tpjjy ^> t c , adjust quasi-statically to the changing 
FUV flux. Long-lived, small, pressure-confined clumps advecting into FUV-irradiated regions 



-13- 



provide a prime example of this slow evolution. Such clumps develop a warm PDR surface 
which slowly gets warmer, thicker and less dense. The clump and the heated PDR surface 
evolve in near pressure equilibrium, with the density of the surface decreasing with the rise 
in its temperature, such that npp>pTpp,p = n c T c = pressure of the ICM (assuming thermal 
pressure dominates on these small scales). At any given time, the evolutionary state of a 
slowly heated clump is the same as the final state of an impulsively heated clump exposed 
to the same local FUV flux. However, a gradually heated clump constantly adjusts its 
PDR surface to maintain pressure equilibrium as the PDR slowly heats. In contrast, an 
impulsively heated clump undergoes shocks and supersonic expanding flows in its attempt 
to attain equilibrium. (We note that there is no final equilibrium state if the interclump 
pressure is zero). Rapid photoevaporative mass loss can cause an acceleration due to a 
rocket effect in very small clumps and move them back into the shielded cloud. Pressure- 
confined small clumps, on the other hand, are not rocketed as they evolve quasi-statically. In 
a given FUV field and at a given ICM pressure, the final equilibrium state of a small clump 
is either a completely photodissociated, warm (Tpdr), and expanded (but with P = Picm) 
region; or it may consist of a shielded cold molecular core with the same temperature, density 
and pressure as the original clump, but with a heated and more diffuse protective surface 
layer of PDR material. In §5.3 we provide analytic solutions for these equilibrium structures. 
However, in §5 and §7 we primarily focus on the interesting physical processes which occur 
for clumps that are impulsively heated. 



5. Analytical models for FUV photoevaporating clumps 

We first provide approximate analytic solutions for the time evolution of a clump sud- 
denly exposed to FUV radiation, using a few simplifying assumptions. As described in §2, the 
nature of the clumps themselves is not well understood, and we consider two simple analytic 
models that qualitatively correspond to the two main interpretations. Turbulent unbound 
clumps are modelled as constant density structures in vacuum. A clump in vacuum disperses 
on a sound crossing timescale, in the absence of FUV illumination. Pressure-confined clumps 
are also treated as being of constant density, but have a surrounding ICM with equal pres- 
sure. In our models, the ICM only serves to confine the clumps and does not itself evolve or 
otherwise affect the clump evolution. We seek solutions for the time evolution of the mass, 
size and mass loss rate of a photoevaporating clump. 
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5.1. Initial conditions and assumptions 

Clumps are assumed to be dense, small spheres of gas, supported by thermal, turbulent 
and magnetic pressures. We assume that the magnetic field B scales with a constant power 
of the density, so that the magnetic pressure is 

P B = B 2 /8n oc n 1 . (2) 

The initial ratio of turbulent and magnetic pressures, Pjvt and Pb, are set by two dimen- 
sionless parameters, 

a = P N t/Pt, 13 = Pb/Pt (3) 

where Pt is the initial thermal gas pressure in the clump. A clump may have significant 
turbulent support and observations of clumps in star-forming regions (Jijina et al. 1999) 
indicate values of a ranging from ~ in cold, dark clouds to > 2 in regions of massive 
star formation. There is observational evidence to suggest that turbulent support in clumps 
decreases on the smaller scales of star-forming clumps (r c0 < 0.1 pc), where a < 1 (Goodman 
et al. 1998). Measurements of magnetic fields and hence (3 and 7 are difficult, but present 
observational data suggest a wide range of values for f3, from ~ to a few (Crutcher 1999). 
The value of 7 (~ 1 — 2) describes how the B field responds to a relatively sudden change 
in gas density. For a frozen, uniform magnetic field threading a spherical clump, B oc n 2 / 3 
or 7 = 4/3. For our standard case we adopt values of a — 1,(3 — 1, and 7 = 4/3. We later 
discuss the sensitivity of our results to these choices for a, (3 and 7. As discussed in §4, 
for clumps confined by an ICM and advecting into FUV-illuminated interclump zones, the 
assumption of tpuv/tc < 1 is equivalent to r c0 > (0.3 — 1)X PDR . Therefore, in this case we 
treat only clumps above a certain minimum size, which is of order the size of the thickness 
of the heated PDR of the ICM. We do not explicitly include self-gravity in our analytical 
model, and only qualitatively estimate the effects of gravity on clump evolution. We treat 
gravity more quantitatively in the numerical models. 

The FUV radiation field as seen by the clump is spherically asymmetric, and stronger 
at the surface facing the source. For simplicity, however, we assume a spherically symmetric 
FUV field when calculating the internal dynamical evolution and the thermal structure of the 
evaporating clump. Typical albedos of interstellar dust are ~ 0.5 and some of this scattered 
radiation is directed backwards (dependent on the phase factor g, see Henyey & Greenstein 
1941), resulting in a backscattered FUV flux on the shielded clump surface about 0.1 — 0.2 
times that on the directly illuminated surface (e.g., Desert et al. 1990, Hurwitz et al. 1991). 
Although the posterior intensity is diminished compared to that on the front surface, the 
PDR surface sound speed is not a sensitive function of the incident FUV flux G (Kaufman 
et al. 1999). For field strengths G ~ 10 3 ~ 5 near typical H II regions and densities n > 10 3 ~ 6 
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cm -3 characteristic of clumps there, even a back flux 10% of that on the clump surface facing 
the source would heat the gas to nearly the same sound speed (cpdr = 3 km s" 1 on the 
source side, cpdr = 2 km s _1 on the shadowed side). Therefore, the assumption of spherical 
symmetry is as good as ignoring the flux of FUV photons on the backside when calculating 
the thermal structure and internal dynamics. However, the asymmetry of the flows from the 
front and back surface will cause a rocket effect on the clumps and induce motion away from 
the source, as first recognized by Oort & Spitzer (1955). We discuss momentum transfer 
due to mass loss of the clumps in §7, and estimate sizes of clumps that gain velocities high 
enough to move out of the PDR and back into the molecular cloud. 

In our analytical models (but not our numerical models) we assume that the transition 
from the warm outer PDR layer to the cold molecular core of a clump is very thin, so that 
the clump can be modelled as an isothermal core with sound speed c c , surrounded by a 
PDR envelope of sound speed cpdr- The gas in the core always remains isothermal in our 
analysis, even when shock-compressed to high densities. For our pressure-confined clumps, 
it is assumed that the ICM is of low density and does not absorb any FUV photons. In 
order to obtain analytic solutions, we do not fully treat transient phenomena and emphasise 
quasi-steady state aspects of the flow. 



5.2. Analytic model for unconfined clumps 

We begin by considering the evolution of a simple configuration: a cloud consisting of 
clumps that are spheres of radius r c o and constant initial density n c o, immersed in a vacuum. 
We refer to such clumps hereafter as "turbulent clumps" . These are dynamic structures that 
in the absence of FUV heating expand at roughly their sound speed, c c . When such a clump is 
exposed to an FUV field, the photon flux is attenuated by dust and instantaneously heats an 
outer column of N ~ 2 x 10 21 cm -2 (thickness So = N /n c o) to relatively high temperatures 
and sound speeds (for example, near H II regions with G ~ 10 4 ~ 5 and n c0 ~ 10 5 cm" 3 , 
T ~ 1000K and c PDR ~ 3 km s" 1 ). 

The evolution of the clump is determined by its initial radial column density from the 
centre outwards, n c or c0 , and the ratio v of sound speeds in the FUV- heated region and the 
cold clump material, 

v = c PDR / c c (4) 

The parameter v can be thought of as a measure of the strength of the FUV field incident 
upon the cold clump. We introduce the dimensionless parameter i] c0 , 



Vco = n c0 r c0 /N = r c0 /S Q , 



(5) 
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for the initial column density to the centre of the clump. Note that 7] c0 is a measure of the 
mean column density through a clump. 

For turbulent clumps with r) c0 < 1, the initial clump column density is less than N , and 
the entire clump is immediately heated and photodissociated by the FUV flux. Effectively, 
therefore, the FUV flux accelerates the expansion of the cold clump by heating it throughout, 
decreasing the expansion timescale by a factor given by the ratio of sound speeds at the two 
temperatures, c c /cpdr, or 1/V. For clumps with t] c q > 1, the FUV flux heats an outer shell of 
gas to a higher temperature, thereby increasing its pressure. If the pressure of the outer PDR 
shell is sufficiently high, it drives a shock to the centre of the clump rapidly compressing it, 
and the compressed clump proceeds to evaporate on a somewhat longer timescale than the 
expansion timescale in the absence of an FUV field. On the other hand, if the pressure in the 
heated outer PDR shell is low, the shock does not make it to the centre of the clump. The 
shock dissipates, followed by an expansion of the clump, until finally the photoevaporative 
flow halts the expansion and proceeds to shrink and compress the gas. There are thus two 
distinct evolutionary scenarios for photoevaporating turbulent clumps with i] c0 > 1. 

For brevity in the main text, the details of our analysis and the derivations of the results 
are presented in the appendices to the paper. We now discuss our solutions and describe the 
physics of FUV-heated clump evolution and refer the interested reader to the appendices for 
the full analysis. 

Clumps exposed to a given FUV field, measured by the parameter u, are either com- 
pressed by shocks or expand to adjust their column density to a critical value, r] crit . This 
critical column density, r] cr it can be derived from conservation of mass and from the condition 
of pressure equilibrium at the surface of the cold clump gas (see Appendix A for details), 
and is defined by the relation 



Recall that a and (3 are the initial ratios of turbulent to thermal pressure and magnetic 
to thermal pressure, respectively, and that the magnetic pressure scales as n 7 . For typical 
values of these parameters, the second term on the LHS is negligible and equation (6) can 
be used to define a photoevaporation parameter, A, as 



For clumps with initial column densities 7] c0 = rj cr i t , the photoevaporation parameter A = 1. 
Figure 2 shows an rj — v parameter plot with equation (6) depicted for the standard case 
with a — l,/3 = 1, and 7 = 4/3, and for two possible extremes of these parameters. The 





A 



(I + a + (3)(r] c0 - 1) 
2(2z/ 2 + a) 



(7) 
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result is not very sensitive to the choice of the parameter 7, and r] crit only depends on the 
sum of a and (5. For the standard case and with r] c0 ^> 1 and v >> 1 which is typical of many 
PDRs, equation (6) can be approximated as r] crit ps 4z/ 2 /3. The evolution of FUV-heated 
clumps is determined by where they initially lie in the 77 — v parameter space. 
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Fig. 2. — Logarithmic parameter plot of the ratio i] c q of the initial column density through 
the clump to the fiducial column iVo, and the ratio v of the sound speed in the heated 
PDR gas to that in the cold clump gas, for impulsively heated clumps. The critical column 
density r\ crit is plotted against v for different values of a and f3, with 7 = 4/3. Also shown 
are representative trajectories of column density evolution for clumps with initial column 
densities less than, and greater than, r] crit . Note that, in this case, the abscissa refers to the 
instantaneous column density through the clump. At t — 0, v — 1 and the clump has a 
column i] c q. With the turn-on of the FUV field, v increases rapidly. A clump with rj c0 < r] crit 
is then compressed by a shock to i] crit and moves right on this parameter plot and then 
proceeds to photoevaporate steadily. A clump with t] c q > i] cr i t expands so that its column 
density becomes equal to r] cr it and continues to photoevaporate. The initial turn-on of the 
FUV field (dotted line) moves the column slightly to the left because a column N is quickly 
dissociated and lost from the clump. Gravity is ignored in this figure (see Figure 5). 
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Clumps with an initial column density r] c0 < r] crit evolve through a shock-compression 
phase, until their dimensionless column density increases to r] crit (see Appendix B). Clumps 
with higher initial column densities expand until their column density decreases to r] CTit (see 
Appendix C). The evolution of clumps in two representative cases is qualitatively shown in 
Figure 2 through their trajectories in 77 — v space. Once the clump column density attains its 
critical value r] crit , further evolution is similar in both cases, and the clump column density 
at any later time remains constant at this critical value. We note that this would imply 
that observed clumps in a given PDR should typically have the same column density, the 
critical value r] crit given by the ambient FUV field. This is valid for clumps photoevaporating 
into a vacuum and only violated in the transient early stages which are typically very short, 
of order t ~ t c q/cpdr- From equation (6), clumps at 10K heated on their surfaces to 100 
K (or v ~ 3 ) should have columns r] cr i t No ~ 3 x 10 22 cm~ 2 , and for T ~ 1000JT or 
v ~ 10,?7 cr j t iVo ~ 3 x 10 23 cm -2 . As photoevaporation proceeds, and clumps lose mass, they 
shrink in size and their average density increases, so that the column through the clump 
remains the same, n c (t)r c (t) = r] crit No. 

Shock compression In this case, the initial column density is larger than Nq, but smaller 
than i] crit NQ. The FUV can only penetrate through an outer shell of thickness 5q = Ao/n c o < 
r c0 . This shell is heated to a higher temperature and the corresponding increase in pressure 
causes the outer shell to expand at cpdr radially outwards at the edge and inwards into 
the clump at the inner boundary. Figure 3 shows a schematic diagram of the evolution of a 
shock-compressed clump. 
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Fig. 3. — Schematic diagram of evolution of shock-compressed clump. The cold clump is 
indicated by the lightly shaded region, the dotted region represents the heated warm PDR 
surface of the clump, and the dark shaded region shows shock-compressed gas. Panel c shows 
the propagation of the shock front into the clump and panel d depicts the evolution of the 
clump after being compressed by the shock into a small dense core. 
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As the shell expands, the column density N through it is maintained by further pen- 
etration of FUV photons, and heating of more cold clump material. The pressure at the 
base of the PDR shell is initially higher than the cold clump gas. This pressure difference 
drives a shock wave into the clump, which rapidly propagates to the centre. Behind the 
shock the cold clump gas is compressed to a pressure approximately equal to the pressure at 
the base of the PDR shell. The shock travels inward at roughly cpdr and the entire clump 
gets compressed in a time t s ~ r c (0) / cpdr, where r c (0) = r c0 — 5 is the initial radius of the 
cold clump when the FUV flux is turned on (Fig. 3b). The shock reaches the centre and 
compresses the clump to a radius r c (t s ) = r s . The shock-compressed gas forms a high-density 
core, such that the pressures in the core and at the base of the PDR flow are approximately 
equal. Further mass loss from the clump is from the surface of this core. The time evolution 
of the mass and radius of the clump for t > t s are given by the following equations for our 
standard case, where a — 1, (5 — 1, and 7 = 4/3. For the general solution, refer to Appendix 
B. 



10r/ c0 z/ 9/4 ft 




dm c f m c0 \ f r c (t)\ 61s 



dt \ t c J \ r c0 J r] c0 



(9) 



m c (t > t.) = m c (0) (^y) 9/4 (10) 

Here q = (A/3) 1 / 3 , r s « gr c (0)(l - 3(1 + q)/(r] c0 - I)) 4 / 9 and t > t s = r c (0)/c PDR . The 
PDR/core interface propagates inward as the warm gas evaporates and, finally, the entire core 
is transformed to heated, photodissociated, expanding PDR material. Using equations (8) 
and (10), Figure 4 shows the clump radius and mass as a function of time for a case where 
the initial column through the clump is given by t] c q = 10, and v — 5 represents the incident 
FUV flux. This is a typical column density for clumps in molecular clouds, (see §2) and 
for a clump of mass m c , corresponds to a density n c0 = 3 x 10 6 (M Q /m c ) 1//2 cm" 3 . In this 
example of v — 5, the FUV field might heat cold clump gas with initial temperature of ~30 
K to about ~750 K, or ~ 10K gas to about 250 K. 
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Fig. 4. — Time evolution of radius (panel a) and mass (panel b) of cold clump gas for a 
clump with i] c0 — 10, v — 5 that undergoes shock compression. The mass and radius are 
given as ratios to the initial mass and radius, and time is in units of sound crossing times 
in the initial cold clump. The solid line is the analytical result, showing an initial rapid 
decrease in clump radius followed by a slower phase of clump radius evolution. The dashed 
line is the result from the 1-D hydrocode and is seen to agree well with the analytical curve. 
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Also shown in Figure 4 are the results of a more detailed hydrodynamics code simulation 
(described in the next section) for the same parameters. There is very good agreement 
between the two results, in spite of the various simplifying assumptions made in arriving 
at the analytical solutions. The numerical solution for r c {t) oscillates with time due to the 
dynamics of the shock compression which we have neglected in the above analysis. Our 
analytic solution assumes a quasi-steady state, whereas the numerical code shows the effect 
of overshooting the steady state solution and rebounding. 

The clump evolution goes through two distinct phases, which is apparent from the plot 
of clump radius with time (Fig. [4a]). The clump initially shrinks rapidly with time, as the 
shock compresses the clump in a time t s . The average velocity vj, with which the radius 
decreases is 0.65cpp,p for the chosen parameters. After the entire clump is compressed by 
the shock, it decreases in size more slowly, now due entirely to the mass loss from its surface. 
As discussed earlier, the column density in the compressed clump remains constant, and the 
clump gets denser as its size decreases. In Figure 4b, the mass of the clump is also seen to 
decrease with time as the entire clump gets photoevaporated in about a sound crossing time; 
again, the analytical results closely follow the result of the numerical hydrodynamical code. 
The cold clump mass at t = is less than the initial clump mass because it excludes the 
mass contained in the outer column N which is instantaneously heated by the FUV flux to 
the PDR temperature. 

Using the results of the above analysis, we can obtain simple estimates of photoevapo- 
ration timescales for a clump exposed to incident FUV radiation. We define the photoevap- 
oration timescale, tpE, as the time for the radius of the clump to shrink to zero. The time 
tpE, can be easily determined from setting the LHS of equation (8) to zero, 

t PE « Otrffiv-VHc (11) 

or 

*- Klo (iow) fcj \-^—> b^rJ years (12) 

Larger clumps are longer-lived compared to smaller clumps in a given PDR environment. 
Photoevaporation timescales of turbulent clumps in PDRs in typical star-forming regions 
are thus of the order of a few clump sound crossing times (see §7). Paradoxically, the sudden 
turn-on of FUV radiation on a clump in a vacuum can increase its lifetime, even though the 
FUV heats the surface and increases the flow speed from the surface! For turbulent clumps 
with lifetimes ~ t c in the absence of an FUV field, exposure to the FUV field results in 
a "pressure-confined" period where the clump is compressed due to pressure of the heated 
surface PDR layer. This compression reduces the area of the photoevaporating surface and 
can extend the lifetimes of these clumps to a few t c . 
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Collapse of clumps driven by shock compression The FUV driven shock wave that 
compresses a clump to very high densities, may render it gravitationally unstable to collapse. 
Thus star formation may be triggered in a previously stable clump. A simple estimate of 
parameters leading to clump collapse can be made without explicitly including self-gravity in 
the equations. We use our solutions for the radius of the shock-compressed core (Appendix 
B) and compare this with its Jeans length, rj. If r s < rj, the clump is stable, otherwise 
it undergoes gravitational collapse. It should be noted here that a clump supported by 
magnetic pressure with 7 > 4/3 will always be stable to collapse regardless of the external 
pressure, as the magnetic pressure increases faster during compression than the gravitational 
energy density. Our discussion here is thus applicable for clumps with insignificant magnetic 
fields, or with a magnetic equation of state where 7 < 4/3. We now solve for the collapse 
criterion as a function of initial clump parameters and the FUV field strength as measured by 
v, for an illustrative case where a — 1, (5 — 1, and 7 = 1. The radius of the shock-compressed 
core is given by 

—Hwt^T (i3) 

The Jeans length of the compressed core is given by 



For collapse r s > rj which yields 

-> 1 + i(£Sf) ( " 2 + irl/3 (15) 

The factor n c0 c 2 in equation (15) is the initial thermal pressure in the clump, and for clumps 
in PDRs is typically of the order of 10 6-7 cm -3 K. For n c0 c 2 = 10 6 cm~ 3 K, equation (15) 
becomes 

49 

Vco>l+ p - 1)1/3 (16) 

Figure 5 shows the collapse criterion(Eq.[16]) on the rj — v parameter plot. Clumps with 
initial column densities greater than ~ 40A r and with initial thermal pressures of 10 6 cm~ 3 
K cannot support themselves against gravity and are unstable to collapse even in the absence 
of an FUV field. This region of instability is to the right of the vertical dotted line in the 
plot. Clumps to the left of this line are initially stable, but sufficiently strong FUV fields 
may drive shocks that compress them and trigger collapse. The FUV fields near OB stars 
may thus trigger star formation in previously stable clumps, and increase the star formation 
rate in clouds, in cases where 7 < 4/3. 
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Fig. 5. — This figure shows the r] — v parameter space and demarcates the regions where 
gravitational collapse can be triggered by shock compression (for 7=1). The dashed line 
indicates the column densities of clumps, with thermal pressure 10 6 cm" 3 K, which are 
initially unstable to collapse. Clumps with this pressure and columns to the left of the line 
are initially stable. The regions to the right of the solid line are where clumps are driven to 
collapse by shocks induced by photoevaporation, for initial thermal pressures in the clump 
Pt = 10 6 . For higher pressures P c0 , both the dotted (rj oc -P^ 2 ) and solid lines (77 oc 
move to the right. 
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Clumps with an initial expansion phase Clumps with initial column densities r} c0 > 
r] crit as given by equation (6) develop very thin (5 <C r c0 ) PDR shells when they are first 
exposed to FUV radiation. This very thin shell is initially at a high pressure and drives 
a shock into the clump on the inside edge of the shell while the outside edge of the shell 
expands at roughly cpdr just as is the case for shock- compressed clumps. However, the 
shock stalls before it reaches the clump centre as the PDR pressure rapidly declines due to 
the shell expansion. The PDR pressure and clump pressure become equal before the shell 
thickness is comparable to the clump radius. We summarize below the detailed evolution of 
initially expanding clumps; details appear in Appendix C. 

Because of the high PDR sound speed and relatively small thickness, the PDR pressure 
initially drops faster than the pressure in the cold clump. The pressure in the cold clump 
eventually becomes higher than that in the surrounding expanding PDR shell, and the cold 
clump gas expands at its sound speed, c c (as it would in the absence of FUV radiation 
since the interclump medium is assumed here to be a vacuum). This expansion takes place 
until a time t e at which the density (and pressure) in the clump gas drops to that of the 
outer PDR layer. For t > t e , FUV photons begin penetrating into the formerly shielded 
cold clump gas. The warm PDR gas now confines the clump gas, and further expansion of 
the cold clump is halted. The cold clump proceeds to lose mass gradually and shrinks due 
to photoevaporation, with an evolution similar to the final solution for shock-compressed 
clumps. The radius and mass of the clump at times t > t e and for a — 1, (3 < 1, 7 < 4/3 are 
given by 

r c (t>t e )=r c (t e )-^^(t-t e ) (17) 

where r c (t e ) = r c (0) + c c t e , 

-^-^)((l^f->) 

and 

mS>te)=m J^±ll(lMY (19) 
Vco V. r c0 J 

Figure 6 shows the change in radius and mass of a clump with initial parameters r) c0 = 
100, v = 3. The clump expands out initially and loses a small fraction of its mass during 
this phase. After its density has dropped, so that there is pressure equilibrium between the 
cold clump gas and the heated layer immediately surrounding it, expansion is halted. As 
photoevaporation causes mass loss off its surface, the clump shrinks in size, and is completely 
heated and photodissociated in about 4 crossing timescales. The numerical results are also 
overlaid on the plot, and there is reasonable agreement. The photoevaporation timescale for 
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initially expanding clumps can be denned analogously to shock-compressed clumps as (see 
Appendix C, Eq.(C7)) 

^-^-K^)(") 2 te)- rs (20) 
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Fig. 6. — Time evolution of radius and mass for an initially expanding clump, with 7] c0 = 
100, v — 3. Radius and mass are given as ratios of their initial values and time is in sound 
crossing time units, t c . The figure shows the analytical result (solid line) and that from the 
numerical code (dashed line). The clump radius increases in the beginning as the clump 
expands and then shrinks as the clump loses mass due to photoevaporation. 
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Therefore in this regime of low is, the photoevaporation timescale tpE is of order several 
sound crossing timescales, and is proportional to v. 

5.3. Analytical model for clumps pressure-confined by an interclump medium 

For a cloud model in which clumps are structures in pressure equilibrium with an ICM, 
the evolution of FUV- heated clumps is somewhat different. The heated shell or PDR can 
now no longer expand indefinitely, but is eventually confined by the pressure of the ICM. 
The evolution of the clump is again determined by i] c q and v. Clumps with an initial column 
density r] c0 < N , are completely heated and photodissociated instantly. Such clumps expand 
until the PDR gas reaches the interclump pressure. Clumps with higher initial column 
densities again are only heated on their surfaces to a column Nq, and their further evolution 
depends on the turn-on time of the FUV field tFuv, relative to t c . 

If the clumps are heated impulsively, tpuv ^ t c , pressure-confined clumps initially 
evolve similarly to the turbulent clumps discussed earlier. Clumps with i] c q < T] cr u are thus 
shock-compressed, and mass flows are set up which photoevaporate the clump. However, the 
outflow runs up against the interclump pressure and eventually reaches pressure equilibrium 
with its surroundings. During this time, the clump may either be completely heated through 
and get transformed into a sphere of warm PDR gas, or may only be partially photoevap- 
orated. The final configuration for partially photoevaporated clumps consists of a small 
remnant cold clump, and an extended warm PDR envelope surrounding it and protecting 
it from the FUV flux. The clump, PDR envelope and interclump medium are in pressure 
equilibrium. Clumps with i] c q > i] crit are not much affected by the FUV field. There is no 
expansion of the clump due to the presence of the confining ICM, and the initially heated 
thin PDR shell expands to form an expanded, but still thin, protective PDR layer. Con- 
tainment of the PDR shell by the ICM prevents further penetration of FUV photons into 
the cold clump and there is no additional heating. For these clumps, photoevaporation is 
relatively unimportant and they retain large fractions of their initial masses. 

Clumps that are heated with tFuv > t c evolve quasi-statically, steadily adjusting them- 
selves to the pressure of the heated PDR shell. The evolution of these clumps can be deter- 
mined from simple steady-state equilibrium considerations. In their final configuration, the 
clumps may retain a fraction of their initial cold gas, surrounded by a warm PDR envelope, 
in pressure equilibrium with the ICM. We consider two extremes in density profiles of the 
clump, a constant density clump and a truncated isothermal sphere (n oc 1/r 2 ) in pressure 
equilibrium with the ICM. 
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First, we solve for the evolution of a clump of constant density. From conservation of 
mass, 

n c0 r c0 3 = n f r f 3 + n PDR (r PDR -r 3 f ) (21) 

where rif and updr are the number densities in the final cold clump core and the PDR, and 
r c o,r/, and rpp> R are the radii of the initial clump, remnant clump, and the PDR envelope 
respectively. As the final configuration is a remnant clump in pressure equilibrium with the 
PDR and the ICM, the final density of the core rif = n c0 , and 

n c0 c c 2 +an c0 c c 2 +(3n c0 c c 2 = np DR cp DR 2 +anp DR c 2 +f3n c oc 2 {n PDR /n c oy (= Pressure in the ICM) 

(22) 

Further, we know that the column density through the PDR envelope is N , as the remnant 
clump is shielded from heating. Therefore, 

np D R,{rpDR - rf) = N (23) 

Using equations (21) and (23) to eliminate updr and rpp>R, equation (22) can be written as 

^ = {l + a + ^- K l ^-a (24) 

where 

e=(^-^J-^) 1/2 -^/ (25) 

and rjf = n c0 rf/N . For complete photodissociation rjf = 0, or £ = rf c ^ . For the standard 
case a — 1, f3 — 1 and 7 = 4/3, and for ?7 c0 > 1, v > 1, equation (24) can be used to define 
a critical column density for complete photoevaporation. 

7] cr ~ 0.48z/ 4/3 (26) 

All clumps with an initial column density lower than i] cr are thus eventually completely 
photodissociated. Alternately, a clump with an initial column density r] cr has to be exposed 
to an FUV field (v) greater than or equal to that given by equation (26) to be completely 
photodissociated. 
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Fig. 7. — The evolution of column density 77/ of the clump against the parameter v. As soon 
as the clump is exposed to an FUV field, an outer column N is instantaneously heated, and 
thus for v marginally greater than 1, at t — 0, the final column density rjf = i] c0 — 1. The 
figure shows how the final column density of the equilibrium clump configuration decreases 
with an increasing FUV field. Also shown are the contours depicting the percentage of mass 
lost from the clump as it photoevaporates. 
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Figure 7 shows a plot of the final column density of a clump, 77/, as it is exposed to 
various FUV fields, measured by is, for a — 1, (3 — 1 and 7 = 4/3. The figure shows 
trajectories followed by the clump column density as the FUV field is increased, as might 
happen when the clump slowly emerges into the PDR region from the shielded interior of 
the molecular cloud. Note that the abscissa is rjf and not the initial column density r] C Q. 
The outer column of cold clump gas N that is initially heated is thereby discounted, and 
as soon as the FUV field is "turned-on" at time t — 0, rjf — rj c0 — I. As the FUV field (v) 
increases, the clump begins to lose mass due to heating of the outer shell, and its column 
density begins to decrease with increasing v. As v approaches the critical value for complete 
photodissociation for that particular i] c0 , the column density rapidly decreases with v and 
the clump photoevaporates completely. If the final strength of the FUV field is less than 
that required for complete photoevaporation, the clump column density remains fixed at a 
point on its trajectory, as determined by the value of the parameter v. Therefore, given 
a clump's initial column density and the strength of the local FUV field, the final column 
density of the clump can be determined. Equivalently, the observed rjf and v can be used 
to determine the initial column i] c0 of the clump. The total mass loss of the clump during 
its evolution can also be evaluated from equation (24). Figure 7 also shows the contours for 
the mass loss being 5%, 50% and 95% of the initial mass, for different values of rjf and v. 

The evolution of a clump with an isothermal density profile can be determined analo- 
gously. The initial density distribution in the clump is given by n(r) = n cs0 r^ /r 2 , where 
n cs0 is the density at the clump surface, and the initial mass of the clump is 47im H n cs orl . 
We assume that the density in the heated PDR shell in the final configuration of the clump 
is constant. Utilizing equations( 21-23) for this case, we obtain an equation very analogous 
to equation (24) where now 



Vcso — n cs or c0 /N and rjf — n cs orf/N . An effective average column density through the 
clump can be defined by the ratio of the mass to area of the clump, f) c0 = m c o/(vrr^ m^A^ ) = 
4?7 cs o. The evolution of a clump with an isothermal density profile is thus qualitatively similar 
to that of a constant density clump, but with a different effective column density for the 
same mass. These solutions apply to clumps that are massive enough for gravity to be 
important in determining their density structure, and therefore have appreciable central 
density concentrations. 




(27) 
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6. Numerical code and results 

In the above analysis, simplifying approximations were made, which are now relaxed in 
a numerical code to obtain complete and more exact solutions (see Appendix D for details). 
One of the major differences between the numerical code and most of the analytic equations 
is the inclusion of gravity in the numerical code. 

6.1. Impulsively heated clump in a vacuum 

We first discuss the results of a numerical simulation of a clump generated by turbulence 
and heated impulsively. This case is hereafter referred to as Case V (for "vacuum"). A 
vacuum boundary condition is used so that the clump is free to expand, even in the absence 
of any heating. Initially, the thermal, turbulent and magnetic pressures in the clump are all 
assumed to be equal (a — (5 — 1), and the magnetic pressure scales with the density to the 
power 7 = 4/3. The FUV field is turned on instantly (t FUV =0). The clump is assumed to 
be a constant density sphere, with a column density given by i] cQ = 10, and in a FUV field 
characterized by v = 10 in the outer layers. The value of i] c0 corresponds to that typical of 
clumps in clouds, for example, m c = 0.8M Q , n c = 2 x 10 5 cm -3 , and r c = 0.03pc. Clumps 
exposed to FUV fields with Go ~ 10 5 are heated to Tpor ~ 1000K on their surfaces; if 
their shielded centres are T c ~ 10JT, then the ratio of sound speeds v ~ 10. These initial 
conditions imply the clump lies in the shock compression regime of Figure 2, and we expect 
from our analytic solutions that it should photoevaporate in less than one sound crossing 
timescale t c (from eq. [11]) 
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Fig. 8. — Density, velocity, pressure and temperature as a function of radius at times t = 
0.05,0.1, and 0.15t c . The density, pressure and temperature are scaled to the initial values 
in the clump and the velocity is in units of the clump thermal speed c c . 
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Figure 8 shows the density, velocity, pressure and temperature profiles in the gas as a 
function of radius at three different instants of time, t = 0.05t c ,0.1t c and 0.15i c . As the 
outer layers heat up and the increased thermal pressure causes them to expand, a shock 
propagates into the clump, compressing the cold inner clump gas. The shock strengthens as 
it progresses to the centre of the cloud, as evidenced by the increasing density and velocity at 
the position of the shock front. The central cold gas gets compressed to very high densities, 
mainly due to the convergence of the radially- moving shocked gas. In Figure 8 we only show 
the evolution prior to the shock reaching the centre of the cloud. After the shock reaches the 
centre, the compressed gas rebounds and the cold inner clump undergoes radial oscillations 
as it settles into pressure equilibrium with the warm, expanding PDR outflow in the outer 
layers. The mean density of the compressed clump gas increases with time, as expected from 
our analytical solutions (eq. [B7]). As the clump loses mass, the FUV penetrates deeper and 
the clump is completely photoevaporated in 0.7t c . Our analytical solution (eq. [11]) predicts 
a photoevaporation timescale of 0.5i c , which agrees very well with the numerical solution. 
The photoevaporation timescale derived for a clump with this column density, r] c o = 10, and 
density, radius and sound speed given by n c = 2 x 10 5 cm~ 3 ,r c = 0.03pc, and c c = 0.3kms _1 
respectively, is 10 4 years. 

6.2. Gravitational collapse 

The compression of a clump by strong shocks due to FUV heating can raise the central 
densities by large factors. This, however, did not render the clump gravitationally unstable to 
collapse in the above case because 7 was chosen to be 4/3, which raises the magnetic pressure 
in dense gas sufficiently to prevent gravitational collapse (Chandrasekhar 1961). The scaling 
index of the magnetic pressure with density in clumps depends on the orientation of the 
magnetic field in the clump, and on whether the field is frozen into the neutral gas. For a 
frozen, uniform, unidirectional magnetic field, Pb oc n 4 / 3 . The magnetic field configuration 
in clumps is considerably more complicated (Ward-Thompson et al. 2000), and ambipolar 
diffusion may allow the neutral gas to slip past magnetic flux lines; making it difficult to 
define 7 unambiguously. 

We investigate the possibility of triggering gravitational collapse in clumps through a 
simulation where we adopt 7 = 1, and all other parameters identical to those in Case V. Such 
a clump is expected to undergo gravitational collapse based on our earlier analysis (eq. [16] 
where r] c0 = 10, and Figure 5). 
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Fig. 9. — This figure shows the mass in the clump, core and warm PDR shell as a function 
of time. At collapse, the clump mass drops sharply to zero. Also shown in dotted lines are 
the clump and PDR shell masses for the case where 7 = 4/3. 
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The evolution of the clump is initially similar to that in Case V. The outer surface 
expands, and a shock is driven into the clump centre. As the shock reaches the centre, the 
density of the cold clump gas increases and the clump radius decreases, until the central 
region becomes gravitationally unstable. A collapsing "core" is thus formed. This can be 
seen in Figure 9, which shows the mass of the clump as a function of time. The mass of the 
collapsing core is also shown in the figure, along with the results of Case V. The core mass is 
seen to increase above that contained in the clump just before collapse as the core accretes 
mass from the warm PDR gas. A significant fraction (60%) of the initial clump mass is driven 
gravitationally unstable to collapse, on very short timescales. In the absence of sufficient 
magnetic field support, exposure to strong FUV fields can thus trigger star formation in 
clumps. We also conducted a simulation with a weaker FUV field {y — 5), and found that 
the central regions collapsed in that case too, but the collapsing core contained only about 
10% of the initial clump mass. Although these triggered collapse solutions are interesting 
and instructive, we emphasize that for initially stable clumps with /3 ~ 1, they hold only for 
7 < 4/3, which requires ambipolar diffusion. Since ambipolar diffusion timescales are longer 
than the shock compression timescale, it is likely that realistic clumps require 7 > 4/3 and 
that triggered collapse does not often occur. 



6.3. FUV heating of a pressure-confined clump: t FUV 3> t c 

Pressure confined clumps are modelled as Bonnor-Ebert spheres which are in hydrostatic 
equilibrium. We construct a Bonnor-Ebert sphere by setting up a density distribution for a 
clump with thermal, magnetic and turbulent pressure support and in hydrostatic equilibrium. 
The density in the clump decreases as r~ 2 in the outer regions (similar to an isothermal sphere 
distribution), and flattens out in the central regions. The ratio of central to surface density x, 
determines the stability of the sphere against gravitational collapse, with a value exceeding 
14.4 being "critically unstable" for an isothermal, non-magnetic gas. In these runs with a 
confining ICM, the outer 200 zones are used to represent the ICM. The ICM has the same 
pressure as the surface of the clump and is a constant density medium with a density (and 
temperature) contrast of 1000. As discussed in §4, for clumps confined by an ICM, a slow 
turn-on of the FUV field may be more appropriate. The FUV field in the numerical run is 
turned on in a time tFuv = 5t c , and we increase the temperature on this timescale so that 
the heated outer shell reaches its maximum temperature at 5t c . 

Figure 10 shows the results for a case with varying x and where r] c o = 10, v — 5, 7 = 4/3, 
and a — (3 — 1. For Bonnor-Ebert spheres with x > 6, the outer region with an isothermal 
density profile contains more than half the mass, and the clump evolution resembles that of 
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an isothermal sphere. In the figure, the mass of the clump is shown with time along with 
the analytical results for an isothermal sphere and a constant density sphere with the same 
initial column. The numerical results for three values of the central to surface density ratios, 
x — 1,3 and 6 are indicated. For central concentrations greater than that represented by 
x — 6, the analytical results for an isothermal pressure-confined sphere match the numerical 
results, and the agreement increases as x increases and the Bonnor-Ebert clump configuration 
approaches a 1/r 2 density profile. For x — 1, the clump is a constant density sphere, and the 
corresponding analytical calculations apply. A slight discrepancy can be noted in the figure 
between the two results, and this is due to the differences in the more realistic numerical 
model, i.e. an exponential drop in the PDR shell temperature with column and the inclusion 
of gravity. For intermediate central density concentrations, as represented by the case with 
x — 3, the mass loss rate is in between the two analytical extremes, as expected. 
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Fig. 10. — This figure shows the evolution of mass with time for a slowly heated, pressure- 
confined clump (Bonnor-Ebert sphere) with varying central to surface density ratios, x — 1, 3 
and 6. In this case, i] c0 — 10, v — 5, and tpuv — 5i c . The corresponding analytical solutions 
for a constant density sphere (dotted line) and an isothermal sphere (dashed line) are also 
shown for comparison. 
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7. Acceleration of clumps due to rocket effect 

Spherically asymmetric mass flows from photoevaporating clumps can cause a rocket 
effect on the clump, as first noted by Oort & Spitzer (1955) in their theory on cloud acceler- 
ation by ionizing radiation from OB stars. Clumps exposed to FUV radiation in PDRs see 
a stronger incident flux of photons from the direction of the OB stars, though a significant 
flux due to backscattered radiation from dust on the opposite surface tends to reduce this 
asymmetry (see §5.2). The front and back hemispheres of the clump (with respect to the OB 
star) are thus heated to different temperatures (different v) and the escaping material has 
a higher flow velocity on the star-facing surface of the clump. The net thrust on the clump 
accelerates it away from the star, and is known as the rocket effect. With the help of our 
analytical expressions for the evolution of an impulsively-heated, photoevaporating clump 
in a PDR, we now evaluate a simple criterion which determines the clump size at which the 
rocket effect becomes significant. 

A simple 1-dimensional formulation for the acceleration due to the rocket effect is given 

by 

dv R 1 dm c 

— -^^r ( ^-^ } (28) 

where Vff and Vfb are the flow velocities at the front and back surfaces respectively and 
are assumed equal to the PDR sound speed cpdr at those surfaces. We first compute the 
acceleration due to the flow from each hemisphere separately, take the difference to estimate 
the net acceleration, and then calculate the net velocity v R attained by the clump. The mass 
loss from a hemisphere is given by 27cr c (t) 2 pbCpDR, where p& is the density at the base of the 
flow. Taking only the component of the momentum in the flow along the direction to the 
star, we obtain 

dv R \ C PDR 2 / on x 

= nr c (t) PbCpDR (29) 



dt J y 2 

The net acceleration on the clump is therefore given by 



,2 



dvR = *r fc 
dt m r 



[(pbCp DR )f ~ (pbC PDR )b\ (30) 



where the suffixes / and b denote the quantities for the front and back surfaces of the clump, 
and the clump radius r/ c is determined by the dominant flux on the front surface. As 
the initial expansion/compression phases in the evolution of a photoevaporating clump are 
typically short (see §5), we ignore this initial phase of evolution and assume the clump is 
evolving such that its column density is at its critical value rj crit ~ 4i/ 2 /3, for the standard 
values of a, (3 and 7 (see eq. [6]). Along with the relation for the column through the flow 
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nb(t)r c (t) = 2N for each surface and the mass of the clump m c = 47imHnf c r 3 j- c /3, (see 
appendix, eq. [A2]), we thus obtain for the net acceleration of the clump, 

dv R _ 9c| 
dt ~ 8 r c 

The clump will remain in the PDR of the GMC and be subject to photoevaporation if the 
rocket velocity v R attained in a time tpp, R = X PDR /vj P is less than vjp, or equivalently if 

9 c 2 c 
Sv 2 n 

There is thus a critical size r R of clumps, where for r c < r R they get rocketed back into the 
shielded molecular cloud, and for r c > r R remain in the PDR and get photoevaporated or 
survive to photoevaporate in the H II region. For example, for front and back surface flow 
velocities of 3 km s _1 and 2 km s" 1 respectively, c c = 0.3 km s" 1 and vjf = 0.5 km s _1 
equation (32) gives a critical clump size r R = 0.225X pp>R- Clumps much smaller than this 
will get rocketed back into the molecular cloud as soon as they enter the PDR, partially 
evaporating in the process. We emphasize again that the rocket effect is significant only for 
turbulent clumps which undergo impulsive heating and photoevaporative flows. Pressure- 
confined clumps are gradually heated on their exteriors and do not go through phases of 
evolution with rapid mass outflows, and hence experience an insignificant rocket effect. 

The mass lost by a clump before it gets accelerated out of the PDR can be estimated by 
integrating equation (28) with time, where only the components of velocity and mass flow 
along the direction towards the source are to be considered. Thus the mass of the clump is 
approximately given by 

m c (t) = m c oexp(-2v IF /Acp DR ) (33) 

where Acpp, R is the difference in the flow velocities on the two sides of the clump. For typical 
values of Acpp, R = 1 km s -1 , and vjp = 0.5 km s -1 , in the time it takes for v R to reach vj P) 
the clump loses about 63% of its initial mass. Therefore clumps which are rocketed back 
into the cold molecular cloud interior also lose significant fractions of their mass before they 
attain high enough velocities to keep up with the advancing ionization front. 
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8. Discussion 

We begin our discussion by summarizing our results for both the case where clumps are 
produced by turbulence in regions with insignificant interclump pressure and the case where 
clumps are initially stable structures confined by interclump pressure. 
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Consider first the scenario where clumps are continually generated by turbulence in 
the molecular clouds and the PDR. Initially, the evolution of the clump is determined by 
its column and the strength of the FUV flux. If A > 1, or t) c q > 4z/ 2 /3 for a — 1, f3 — 
1 and 7 = 4/3, there is a very brief period of shock compression but the shock quickly 
dissipates before making it to the clump centre and the clump radius shrinks only slightly. 
Subsequently, the clump expands until A ~ 1, at which point the photoevaporative flow halts 
the expansion and proceeds to shrink and compress the clump gas, maintaining A = 1, or 
a column n c r c = Au 2 N /3. If A < 1, then the shock propagates to the centre of the clump, 
compressing it significantly, followed by an evolution at A ~ 1 identical to that above. 

During the final A = 1 evolution of the shrinking, compressing, evaporating clump, the 
rocket effect may play a significant role depending on the physical size of the clump at this 
time. If r c < t\r, the clumps will be accelerated into the molecular cloud, losing a significant 
fraction of their mass in the process. If tr < r c < Xpdrc c /vif, the clumps will tend to 
evaporate in the PDR region. However, because in the turbulent model clumps are forming 
in the PDR, including the region near the IF, some of these clumps will survive to enter an 
advancing H II region. If r c > Xporc c /vif and an H II region is advancing, then most of 
these large clumps will survive to be engulfed by the H II region. 

The lifetimes of turbulent clumps are of order several sound crossing times of the initial 
size of the clump (Eq. [11] and Eq. [20]). As discussed previously, this is somewhat longer 
than clump lifetimes in the absence of an FUV field. Since turbulent clumps are being 
continuously formed, the steady-state abundance of clumps could actually increase in the 
presence of an FUV field. 

Based on these results, the following predictions can be made for the turbulent scenario 
of photoevaporating clumps in PDRs. 

1. There will be an enhanced population of clumps with columns n c r c m 4u 2 N /3. Note 
that this implies that in comparing this population from one PDR to another, the 
PDR with the lower v (e.g., due to lower Go) will have clumps with smaller column 
densities. 

2. There will be a population of small clumps (r c < r R ) moving with velocity v ~ vjp 
into the molecular cloud at A v > 1. 

3. Compared with clumps in FUV-shielded regions, the clumps in turbulent PDRs will 
be smaller, denser and potentially more numerous. 

4. Similarly, the clumps entering the H II regions will be smaller and denser than the 
molecular cloud clumps, and will have columns peaked at 4i/ 2 iV /3. 
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We note here that in the present analysis, turbulent clumps are geometrically idealized 
as constant density spheres. In reality, turbulently generated structures may also be sheet- 
like or filamentary in nature. Let us consider two simple representations of these geometries, 
disks and cylinders. For these structures to survive immediate photodissociation, the column 
density through their shortest dimension has to be greater than Nq. Disks with larger column 
densities through their midplane establish photoevaporative flows in an FUV radiation field. 
When the outer boundary of the heated gas reaches a distance equal to the radius of the 
disk, the flow diverges and the density in the flow (refer Eq.[Al] and Eq.[A2]) begins to drop 
as r -3 if v oc r. Therefore, the mass loss rate from the disk is essentially similar to that of 
a spherical clump of size equal to the disk radius (Johnstone et al. 1998) Photoevaporative 
flows from a cylindrical structure differ slightly from those off a spherical clump of the same 
radius. The density at the base of the flow is now oc v _1 r _1 oc r~ 2 and therefore, the critical 
column density r) cr i t ~ 2i/ 3 /3. The mass loss rate from the cylindrical clump (following shock 
compression or initial expansion) no longer depends on the radius and is a constant given by 
—2irmHN cpDRl where / is the length of the cylindrical filament, (cf. Eq. [Bll]). In summary, 
for turbulent clumps which are likely to be sheets or filaments, the mass loss from a thin 
sheet or disk of dimension r is best modelled as a clump of radius r. The mass loss rate from 
a disk of radius r and thickness t is similar to that of a sphere of radius r and is oc r, but the 
evaporation timescale is m/(dm/dt) oc r 2 t/r oc rt and hence much smaller than the sphere 
(m / (dm / dt) oc r 2 ) for t <^ r. Filaments of radius r evolve somewhat differently than spheres 
of radius r; their mass loss rates do not decline as r shrinks. However, the evaporation 
timescales evolve similarly for cylinders or spheres since for spheres m/(dm/dt) oc r 3 /r oc r 2 , 
whereas for cylinders m/(dm/dt) oc r 2 l/r°l oc r 2 . 

Consider next the scenario where clumps are pressure-confined structures which are 
bounded by a low-density ICM, and the FUV turn-on timescale is of order the crossing time 
Xpdr/vif for clumps through the PDR. For small clumps with r c < Xpp,pc c /vjF, the FUV 
turn-on time tpuv > ta which means that the FUV quasi-statically heats and expands an 
outer PDR surface on the clump, continually maintaining a pressure Picm in the PDR region 
and in the cold clump region. Again, the evolution of the clump is determined by its column 
and the strength of the FUV flux. If r] c0 < 0.48z/ 4/3 (for a = 1, (3 = 1, and 7 > 1), the 
clump is completely heated and photodissociated as it enters the PDR from the molecular 
cloud. If T] C Q > 0.48z/ 4//3 , the clump shrinks to a smaller cold clump surrounded by a PDR 
shell. Since the evolution is quasi-static, there is insignificant rocket effect on these clumps. 
In this scenario, small cold clumps, protected by PDR shells, can survive passage through 
the PDR and into the H II region. However, the minimum initial column for such a clump 
is r] c0 « 0.48z/ 4 / 3 (for a = 1, f3 = 1, and 7 > 1). 

For large clumps, r c > Xpdrc c /vif, tpuv < t c , which means that shocks compress 
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the clump and strong photoevaporative outflows are initiated as the clump enters the PDR. 
Clumps with 1 < i] c0 < 4z/ 2 /3, are compressed by shocks which propagate to the centre, but 
they finally relax to an analogous situation as above, i.e., a thick PDR shell surrounding 
a small mass cold clump if i] c q > 0.48z/ 4 / 3 and totally photodissociated otherwise. Clumps 
with i] c q > 4z/ 2 /3 have a brief period of shock compression and PDR expansion, but rapidly 
relax to a configuration with a thin PDR shell and most of the initial mass still in the cold 
shielded core. All the large clumps experience a rocket effect. However, as shown earlier, 
large clumps with r c > r R do not experience enough rocket acceleration to prevent them 
from entering the H II region. Therefore, in the ICM confined scenario, a variety of clump 
sizes enter the H II region. However, the small clumps entering the H II region were initially 
clumps of larger size and mass in the molecular cloud, and they enter the H II region with 
large protective PDR shells. 

One can compare the abundance of clumps of various sizes in PDRs with the abundance 
of clumps in the molecular cloud, under the assumption that ICM-confined clumps are long- 
lived (the formation time of clumps exceeds the crossing time Xpdr/vif)- Small cold clumps 
are converted to PDRs but larger cold clumps are converted to small cold clumps with PDR 
shells and this partially replaces the small clump population. Thus, there will be a significant 
drop in the cold clump population with i] c0 < 4z/ 2 /3. However, the population i] c0 > 4z/ 2 /3 
is little affected. Therefore, a variety of clump sizes enter the H II region, but there is a 
suppression of clumps with columns i] c q < 4z/ 2 /3 compared to the distribution in the shielded 
molecular cloud. 

Based on these results, the following predictions can be made for clumps initially con- 
fined by interclump pressure. 

1. Compared with small cold clumps (r] c0 < 4z/ 2 /3) in similar pressure regions of the 
molecular cloud, clumps of the same column density in PDRs will have a smaller 
relative population. Correspondingly, small and large clumps enter the advancing H II 
regions, but the relative number of small clumps is suppressed. 

2. Larger clumps are relatively unaffected, but do experience an initial transient period 
of shock compression and photoevaporative flow. 

3. Small clumps (r] c0 < 4z/ 2 /3) in the FUV-illuminated region will have small masses 
relative to their PDR shells. 

We have shown that small clumps tend to be destroyed before they are overtaken by the 
ionization front of the advancing H II region. This would affect the number, size distribution, 
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and structure of pressure-confined clumps that enter an H II region, and therefore affect the 
evolution of the H II region (see Bertoldi 1989 and Bertoldi & McKee 1990 for a discussion 
of the effect of photoevaporating clumps on the evolution of an H II region). Clumps that 
survive to enter the H II region will thus have undergone significant evolution from their 
initial states. These clumps could be considerably denser if they are still vigorously photo- 
evaporating as they enter the H II region, where they continue to photoevaporate (Bertoldi 
1989, Bertoldi & McKee 1990). We will discuss the propagation of clumps into H II regions 
in a separate paper, with application to the Eagle Nebula. 

We conclude our discussion by using our results to critique the previous inferences of 
dense (~ 10 7 cm -3 ) FUV heated clumps in PDRs and to propose several observational 
consequences of our results. 

Our photoevaporation models provide arguments against the interpretation that dense 
(n ~ 10 6-7 cm~ 3 ) PDR surfaces of clumps give rise to the observed high excitation CO 
lines. The fundamental problem is the propagation of the FUV photons to regions of such 
high density. If densities of ~ 10 6 ~ 7 cm~ 3 are in fact heated to T 3> 100K, then the thermal 
pressures in these regions are much larger than the pressures in the interclump gas. The only 
way this is possible is if: (i) the clumps are gravitationally bound or (ii) the warm surface 
regions are photoevaporating, as discussed in this paper. In either case, the density in the 
clump will smoothly fall off with radius, n(r) oc r~ b where b ~ 2 — 3, until the transition to 
the ICM. However, FUV photons which heat the CO are only able to penetrate through a 
column N ~ 2 x 10 21 cm -2 , where 



This equation reveals the size r of the FUV heated region if a density n(r) is required in this 
region. It implies that for the FUV to heat gas at densities > 10 6 cm -3 , the clumps must be 
very small with size scales of r N q /10 6 cm~ 3 as 2 x 10 15 cm at this density. Gravitationally 
bound clumps have densities n ~ 10 6-7 cm -3 at sizes ^> 10 15 cm (e.g., Shu 1977). In such 
clumps, the FUV would not penetrate to such high densities, but would be absorbed in the 
lower density regions further out. Such small scales as 10 15 cm also argue against unbound 
clumps explaining the high temperature CO emission. We have shown in this paper that 
such a small clump will either be rocketed out of the PDR or photoevaporate in roughly a 
sound crossing time, or ~ 300 years, a time so short as to rule out this possibility. 

The size scale of 100 AU is suggestive of disks around young stars rather than clumps, 
so that one might appeal to photoevaporating disks like the "proplyds" in Orion (e.g., John- 
stone et al. 1998). Here, the enormous mass reservoir of the disk greatly lengthens the 
photoevaporative timescale. However, the area filling factor of these photoevaporating disks 




(34) 
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required to match the CO high-J line intensities in, for example, the Orion Bar region is of 
order ~ 0.1 (Burton et al. 1990). Such a high area filling factor requires a volume density 
of > 3 x 10 5 stars pc -3 , more than an order of magnitude higher than that observed in 
the densest region of the stellar cluster near the Trapezium (McCaughrean et al. 1994). 
Photoevaporating proplyds are thereby also ruled out. 

Therefore, it seems quite unlikely that FUV radiation is heating dense, n ~ 10 6 ~ 7 
cm~ 3 , PDR surfaces of clumps or that photoevaporating proplyds can explain the mid- J CO 
emission. We suspect that the CO emission arises in n ~ 10 5 cm~ 3 surfaces of large (> 10 16 
cm) clumps, and that the chemistry, heating and/or dynamics needs to be modified in PDR 
models in order to better match the observations. 

There are several observational consequences implicit in the scenario of a strongly pho- 
toevaporating (i.e., tpuv/tc < 1) clump. First, the molecular material from the clump 
core is effectively advected out into the PDR region, as has been discussed by Bertoldi & 
Draine (1996), Storzer & Hollenbach (1998), Hollenbach & Tielens (1999). This can have 
strong effects on the chemistry, primarily because H 2 can exist in surface regions where it 
would have been atomic in a stationary, steady-state case. In turn, the higher abundance 
of H 2 can chemically enhance the abundance of other minor species, such as Hg (Bertoldi 
& Draine 1996) and CO and can modify the heating and cooling processes and therefore 
the PDR temperature. Initial results (Storzer & Hollenbach 1998 and Storzer 2000, private 
communication) suggest that advection will enhance the intensities of the mid-J CO lines, 
by creating more warm CO near the surface. Finally, the flow will broaden line widths 
to values Avfwhm ~ Cpdr- However, it should be noted that many common molecules, 
such as CO, HCO+, CS, HCN, CN still exist primarily at column depths iV > N Q from 
the surface, where the advection velocities are low. Broader linewidths should be seen in 
C II(158/xm),0 I(63//m) and possibly the H 2 (2/im) lines. 



9. Applications to observed PDRs in well-studied star forming regions 

We compare our results with available observational data on clump characteristics from 
some well-studied star- forming regions, the Orion Bar, M17SW, NGC 2023 and the PDR 
surrounding the Rosette Nebula. Clump column densities (rj) can be estimated from mea- 
sured or inferred densities and sizes of clumps in these PDRs. Densities are usually inferred 
indirectly from comparing observed line intensities and line ratios with available PDR chem- 
ical models in most of these cases. Clump sizes are ill-determined, as they are usually too 
small to be resolved directly and only upper limits to the sizes are available, with the PDRs 
in Rosette and NGC 2023 being possible exceptions. The FUV field strengths in these PDRs 
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can be estimated more accurately, as the luminosities of the illuminating O/B stars is rea- 
sonably well-determined. We use these field strengths and the dust temperature models of 
Hollenbach, Takahashi & Tielens (1991) to estimate the temperature of dust in the shielded 
clump interior. As the interior gas temperature is closely coupled to the dust temperature, 
we set the two equal and thus determine the sound speed c c in the cold clump gas. The 
temperature to which the FUV field heats the surface clump gas is determined from the PDR 
models of Kaufman et al. (1999), and thus cpdr is obtained, allowing us to estimate the 
value of the parameter v in each region. Below we first determine rj and v for each observed 
region, and then use these parameters, compared to Figures 2 and 7, to infer their nature 
and evolution. 

Orion Bar The Orion Bar has been observed in many atomic and molecular line transitions 
in the sub-millimeter and infrared wavelength regions (e.g. van der Werf et al. 1996, Tauber 
et al. 1994, Young Owl et al. 2000, Marconi et al. 1998). The molecular line emission 
peaks at about 0.03 pc from the IF, and the gas between this dense ridge and the IF is 
mainly neutral atomic hydrogen. The observed line intensities and line ratios appear to be 
best explained by a two-component model for the PDR gas, with small (< 0.02 pc) dense 
(n c ~ 10 6 ~ 7 cm~ 3 ) clumps, required to match the line emission and an interclump medium 
nicM ~ 10 4-5 cm~ 3 which causes the chemical stratification of the edge-on PDR. The local 
FUV field is estimated as G ~ 4 x 10 4 . Such an FUV field incident on clumps with these 
inferred densities would heat their surfaces to temperatures ~ 2000 — 5000 K (Kaufman et al. 
1999). The gas temperature in the interior of clumps exposed to this FUV flux is expected 
to be ~ 50 K (Hollenbach et al. 1991). We thus obtain rj ~ 15 - 150 and v ~ 8 - 10. 

M17 SW Multilevel molecular line observations in CS (Wang et al. 1993, Evans et al. 
1987) and in several fine structure lines (C II, Si II, O I, Meixner et al. 1992) indicate that the 
M17 SW cloud core consists of numerous high density clumps with densities n c ~ 1 — 5 x 10 5 
cm~ 3 . Maps of CO + emission compared with theoretical models suggest that the emission 
is produced in the warm surface layers of PDRs in dense clumps (Storzer et al. 1995) of 
sizes r c ~ 0.1 pc. The local FUV field near the H II region/molecular cloud interface is 
G ~ 8 x 10 4 , which implies PDR temperatures of about 2000K on the clump surfaces and 
cold clump gas temperatures ~ 60K, therefore, v ~ 4 — 9. The column densities of clumps 
are typically rj = 7 — 35, as estimated from typical densities and sizes (Meixner et al. 1992, 
Storzer et al. 1995). 
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NGC 2023 The well-studied PDR in the reflection nebula NGC 2023, at a distance of 
475 pc, is illuminated by a B star and is well-studied (e.g., Wyrowski et al. 2000, Steiman- 
Cameron et al. 1997, Howe et al. 1991). Observations of atomic fine structure lines indicate 
the presence of clumpy gas in the PDR with densities ~ 10 5 cm~ 3 , and an ICM density of 
~ 10 3 cm -3 (Steiman- Cameron et al. 1997). The PDR is situated about 0.1-0.2 pc away 
from the exciting star, and has a width Xpdr = 0.04 pc. The clumps have sizes ~ 0.05 pc 
(Wyrowski et al. 2000), and therefore rj c0 ~ 5. The FUV field strength has been estimated 
at G = 10 4 , and from the PDR models of Kaufman et al. (1999), c PDR = 3km s~ 3 . With 
a clump gas temperature of ~ 20 K for this value of G , we obtain v = 5. It should be 
noted that the level of turbulent energy in the NGC 2023 PDR as estimated from observed 
linewidths is low (Wyrowski et al. 2000), which suggests that clumps are not rapidly re- 
formed by turbulence in the PDR. 

Rosette PDR The Rosette Nebula surrounding the young open cluster NGC 2244, is an 
expanding H II region at a distance of 1.6 kpc. The H II region/molecular cloud interface 
shows a distinct ridge of emission, with substructure on small scales down to 0.5 pc (Schneider 
et al. 1998). The PDR region is very extended, and the FUV photons penetrate deep into 
the cloud with an interclump medium density tiicm ~ 10 cm~ 3 (Blitz 1991, Williams et 
al. 1995). The ICM near the edge of the HII region may even by partially ionized gas as 
suggested by the pervasive Ha emission observed (Block et al. 1992). This region has a 
low FUV field, with G ~ 200 about 15 pc away from the OB star cluster, and dropping to 
lower values of 10-50, about 30 pc away from the cluster (Schneider et al. 1998). Clump 
densities required to match observed CII line intensities are ~ 10 4-5 cm -3 . Cold clump gas 
temperatures are calculated to range from 10 — 15 K. The Rosette PDR thus has t] c q ~ 5 — 35 
and v ~ 5, while further away from the H II region, v ~ 2 — 3. 

Using the above determined values of the parameters i] and v for clumps in the above 
PDRs, we can locate them on the parameter plots of Figures 2 and 7. Figure 11 shows an 
overlay of our model results with inferred data for both the turbulent and pressure-confined 
clump models. The data are consistent with both clump models and, within the limits of the 
uncertainties in arriving at "observed" rj and v, agree reasonably well with our predictions. 
Turbulent clumps are expected to have columns fairly close to the critical value r] crit which 
is ps 4z/ 2 /3. The observed data lie near the rj = r] crit line, indicating that the observed 
clumps may be turbulent in origin and undergoing photoevaporation in the FUV field. The 
exception is NGC 2023. The clumps in this PDR are observed to have negligible turbulence, 
and are thus not expected to have a column equal to i] cr i t . The data are also consistent 
with clumps being pressure-confined. In such a situation, we would not expect to see clumps 
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beyond the 95% mass loss line, as seen in the figure. However, there does appear to be 
a weak trend of increasing rj with u, which if real is not easily explained in this model. 
The clumps in NGC 2023 are probably pressure-confined, and seem to have lost almost 
75% of their original mass through FUV heating. High resolution observations of clumps in 
PDRs, with more accurate determinations of the sizes and densities of clumps would help 
in locating them more precisely on the rj — u plots, and thus be able to distinguish between 
the pressure- confined and turbulent clump models. 
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Log 77 Log 77 

Fig. 11. — Observational data for some PDRs on the rj — v parameter plot for the turbulent 
(panel a) and pressure-confined (panel b) clump models. There is a slight indication that 
the columns lie along the z/ 4//3 lines expected in the turbulent model of clumps, with the 
exception of NGC 2023. 
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10. Conclusion 

We have studied the effects of FUV radiation from young OB stars on the evolution of 
dense structures or clumps in photodissociation regions. We find that the ambient FUV field 
penetrates through the surface of a dense clump and heats this surface layer to high tem- 
peratures, causing mass loss and thereby inducing photoevaporation of the clump. Through 
analytic approximations and numerical hydrodynamical calculations, we determined the evo- 
lution and lifetimes of clumps, subject to various physical conditions and initial parameters. 

The evolution of a clump is mainly determined by three parameters, the ratio of the 
initial column density to the column penetrated by the FUV flux, r] c0 ; the strength of the 
FUV field, which we denote by the ratio of the sound speeds in the heated gas and the 
cold clump material, v\ and the timescale for the turn-on of the FUV field, relative to the 
sound crossing time through the clump, tFuv/t c - We consider the evolution of turbulent, 
impulsively-heated clumps and also pressure-confined clumps in a PDR. Impulsively heated 
clumps with a photoevaporation parameter A ~ 3?7 c o/(4z/ 2 ) < 1 are initially compressed by 
a shock driven by the external high-pressure FUV-heated layer to high densities, and later 
lose mass as the FUV gradually penetrates through the entire clump. Clumps with A > 1 
go through an initial expansion phase after the rapid decay of the shock wave. Both clumps 
with A < 1 and A > 1 evolve toward the A = 1 condition. The final evolution toward 
complete photodissociation and heating occurs at constant column through the clump, the 
density increasing as the clump shrinks. Photoevaporation timescales in all these cases are 
typically a few t c , suggesting clump destruction times of ~ 10 4 ~ 5 years, under typical PDR 
conditions. Slowly heated clumps evolve quasi-statically, and if confined by an external ICM 
pressure, a fraction of the initial cold clump mass may be retained if rj c o > 0A8u 2 . Clumps 
with initial column densities i] c0 < 0A8u 2 are completely photoevaporated by the FUV flux. 

We predict that in the turbulent scenario, observed clumps should all have columns close 
to the critical value for the local FUV field, i] crit ps 4z/ 2 /3. Clumps lifetimes are prolonged by 
photoevaporation and we also expect a higher steady-state abundance of clumps in PDRs as 
compared to the shielded molecular cloud interior. For pressure- confined clumps, there will 
be a decrease in the number of smaller clumps in PDRs, and many clumps will lose substantial 
fractions of their mass to extended warm PDR shells around them. We compared our results 
to observations of some well-studied PDRs, such as the Orion Bar, M17SW, NGC 2023 and 
the Rosette PDR, and find that the data are consistent with both clump models, but perhaps 
favour the turbulent clump interpretation. Clumps entering the H II region around an O /B 
star with a D-type ionization front are expected to have thus undergone significant evolution 
as they pass through the PDR. 
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A. Analytic model: Unconfined clumps with tpuv/t c <^ 1 

We first define several characteristic timescales during evolution. The initial sound 
crossing time through the clump is denoted by t c = r c o/c c . We can similarly define a sound 
crossing time through the warm PDR layer, as t p = Sq/cpdr- Finally, we define the timescale 
for the PDR layer to expand to a thickness comparable to r c0 , t e0 = t c0 /cpdr- 

We relate the density n b at the base of the PDR shell (r = r c ) to N for epochs when 
the PDR shell has expanded to a thickness greater than the clump radius (t > t e0 ). The 
FUV-heated shell of initial thickness 5$ expands as an isothermal flow into vacuum and the 
characteristics of the flow regions are essentially similar to that of the Parker wind solution 
with no gravity (Parker 1958). For simplicity, it is assumed here that once the warm PDR 
region has expanded to 5 > r c or t > t e o, the gas leaves the clump surface at cpdr and the 
velocity of the heated gas increases linearly with radius. Based on our hydrodynamical code 
results, this is a reasonable approximation to within about two clump radii, where most of 
the column density in the PDR gas lies. The clump mass loss rate is determined by the 
penetration depth of the FUV flux, and hence Nq. For a power law density profile in the 
shell (n(r) oc r~ b , b > 2), most of this column arises from the base of the shell, or, more 
precisely, from the region between the base (clump surface) r c and ~ 2r c . From the equation 
of continuity and v oc r, the steady flow density profile of the warm PDR shell is obtained 
as, 

n(r) oc v~ l r~ 2 oc r~ 3 , (Al) 

valid for r c < r < 2r c once the shell has expanded to at least two clump radii (t > t e o). 
Equation (Al) can be easily integrated over r to obtain the column through the warm PDR 
at any given instant of time. Let r c denote the instantaneous cold clump radius and r t the 
total radius of cold clump and warm PDR gas. At t = 0, r c = r c (0) = r c0 — S . Using the 
fact that r t 3> r c for t > t e0 , and as the FUV flux always penetrates through a column A^, 
we have r 

N = f n(r, t)dr w ^n b (t)r c (t) for t > t e0 (A2) 

We next write down the equations of momentum flux or pressure. The pressure in the 
shocked cold clump gas is the sum of thermal, turbulent and magnetic pressures and given 
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by 



P c (t) = n c (t)c 2 c + an c (t)c 2 c + (3n c0 c 2 c (n c (t)/n c0 ) 



7 



(A3) 



where n c (t) is the shocked clump density at time t, n c0 is the initial clump density, a is 
the initial ratio of turbulent and thermal pressures, and f3 is the initial ratio of magnetic to 
thermal pressure in the clump. The pressure at the base of the PDR flow, P b , is given by 



where vji ow is the outward flow at the base. We assume that Vfiow = cpdr, and from 
equations (A2) and (A4), 



The boundary between the "shock compression" and "initial expansion" regions of Fig- 
ure 2 is determined by comparing equations (A3) and (A5). Values of the parameters i] c0 
and v on this boundary are such that the shock dies out just before making it to the cen- 
tre. The pressure in the clump at this instant is about the same as the initial pressure and 
r c (t e o) ~ r c (0). In a time t = t e0 , the outer edge of the PDR shell expands out to 2r c0 , and 
the initial PDR column begins to decrease appreciably due to spherical divergence, allowing 
further FUV flux penetration. 

If Pb(t e0 ) > P c (t = 0), the pressure in the PDR flow is sufficiently high for the shock 
wave to reach the centre of the clump, compressing the clump and raising its central pressure 
so that Pb(ieo) = Pc(teo)- Such clumps lie in the "shock compression" of Figure 2. If the 
initial parameters are such that Pb{t e o) < P c {t = 0), the shock wave rapidly dies out before 
reaching the centre and the clump then stops contracting and begins to expand at a speed 
c c , maintaining pressure equilibrium. These two evolutionary sequences are separated by 
the (r/co, v) boundary obtained setting Pb(t e0 ) = P c (t = 0), or from equations (A3) and (A5), 



For values of r/ c0 and v typical in PDRs such as Orion (r] c0 3>l,z/ 2 3>l,o; = /3 = l,7 = 
4/3) this relation can be approximated as r] c0 ps 4z/ 2 /3. A photoevaporation parameter 
A = 3(r] c0 — l)/2(2z/ 2 + 1) can be defined such that the relation A = 1 demarcates the two 
regions. Clumps with A < 1 are shock compressed and those with A > 1 initially expand. 




(A4) 




(A5) 




(A6) 
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B. Evolution of shock-compressed clumps 

As the flow expands in a shock compressed clump, the clump shrinks by losing mass at 
a rate 

d>7Tl 

— - = -Aixm H n b (t)r c (t) 2 c PDR (Bl) 

Initially, a shock is driven into the clump, but we do not attempt to analytically model the 
propagation of the shock wave in detail. We assume that the shock travels inward at the 
PDR sound speed and the entire clump gets compressed in a time t s = r c (0)/cpoR, where 
r c(0) = T c o — ^o- We first solve for the clump radius, mass and mass loss rate for t < t s , 
during the epoch when the shock propagates to the centre. 

We assume that the clump radius decreases with a constant speed v b (< c P dr) so that 

r c {t) = r c (0) - v b t. (B2) 
The mass of the cold clump decreases with time, and at t s is given by 

m c (t s ) = m c (0) + f (dm c /dt)dt (B3) 
Jo 

where m c (0) = |7rm#n c0 r c 3 (0). From equations (A2), (Bl), (B2) and (B3), we have 

m c (t a ) = m c (0) - 87rm H A r r c (0) 2 (l - v b /(2c PDR )) (B4) 

We solve for v b and thereby calculate r c (t s ). To do this we need to calculate the clump 
compression at t s . The clump gets compressed by the shock to a new density n(t s ) which 
can be determined from the condition of pressure equilibrium at the clump surface. The 
momentum flux conservation equation across the "front" which marks the boundary between 
clump and FUV heated flow is then 

(1 + a)n c (t s )c 2 c + (3n c0 c 2 c (n c (t s )/n c0 y = 2n b (t s )c 2 PDR + an b (t s )c 2 c + (3n c oc 2 c (n b (t s )/n c0 y (B5) 

In moderate to strong FUV fields, the PDR sound speeds and hence the shock velocities 
(assumed equal to cpdr) are high, and the resulting compression is also high. As the magnetic 
pressure Pp., scales with a higher power of the density than the thermal pressure P p , (we 
use 7^4/3) and as the initial magnetic pressure is comparable to the thermal pressure, the 
thermal contribution to the pressure in the compressed gas can be ignored. In the expanded 
flow, the density is low and here the thermal pressure and dynamical pressure dominate. 
Equation (B5) can be simplified to give the density of the compressed gas 

n c (t s ) ((2v 2 + a)n b (t s )\ lh 

(B6) 
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From equations (A2) and (B6), 



no V P ) ^ ■ <B7) 
The mass of the clump at time t s can also be expressed as 

4<7j- 

m c(t s ) = -Yrn H n c (t s )r c (t s ) 3 . (B8) 
From equations (B2), (B7) and (B8) 

m^^ymPcoM^ -J r c (0) 3 ~ 1/7 f 1 - J . (B9) 

At time £ s , the radius of the shock-compressed clump is r s = r c (t s ) = r c (0)(l — v^/cpdb) 
from equation (B2). The unknown parameter is finally determined from equating the two 
expressions for the clump mass, equations (B4) and (B9). The evolution of the clump radius 
is thus known as a function of time for < t < t s , from equations (B2),(B7) and the solution 
for Vb- These general equations are, however, fairly complex, and an approximate solution 
can be obtained for the case a — 1, (5 — 1, and 7 = 4/3 as, 

r s » qr c (0) (l - " (B 1 ) 

V Vco — 1 / 

where q = (A/3) 1 / 3 . 

We now solve for the evolution of the clump for t > t s . After being compressed by the 
shock, the clump continues to lose mass through photoevaporation, with a mass loss rate, 

dun 

= -Aixm H n b (t)r c (t) 2 c PDR = -8nm H N r c (t)c PDR (Bll) 
The mass of the clump at any given time (t > t s ) can also be written as 

From equation (Bll) and differentiating equation (B12) with time, the radius of the clump 
as a function of time is given 



r c (t) 



2 ~ 1/T _ (2 7 -i) 6u Vc0 / ^ c0 -i) \ 1/7 _ fao-ir x ^ 

m) ' (37 - 1) (Vco - l) 2 V 2 (2^ 2 + <*)) \t c ' Vc0 u 



(B13) 
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A photoevaporation timescale tpE for the clump can be defined as the time for the radius of 
the clump to shrink to zero, and from equation (B13) 

tPE = tc ( (*) (h^i) M (^±4) lh + v«zi) . (B14 ) 



r c (0)J V 2 7-l/ Q^Vco \/3(r} c0 -l)J rj c0 u 

For a — 1, (5 — 1, and 7 = 4/3, the above equations determining the evolution of the clump 
can be simplified as 

r s V /4 iO^co^ 9/4 ( t (r) c0 -iy 



'«<«> = ! I^J -3CT^ /4 l"^ll '•«»• < B15 > 

mc(t) = mM (^o)) ' (B16) 

and 

's^-l) 2 / 3(l + g)\ 5/9 , (ifco-1) 



*™ = *cRf U-'^r + ~ ' ■ ( B1 7) 



C. Evolution of clumps that undergo an initial expansion 

Clumps with large initial column densities and in low FUV fields quickly evolve to the 
point where their internal pressures are greater than that of the (expanded) thin heated 
surface layer. They expand out into vacuum (as they would even in the absence of an 
external FUV field) at their sound speed c c for a time t e , until the pressure drops to that 
in the heated outer layer. At t = t e , there is pressure equilibrium and equation (B5) again 
holds. The clump expands to a new radius r c (t e ) = r c (0) + c c t e and during expansion loses 
mass at a rate given by equation (Bl). The mass of the clump at t — t e is obtained as earlier, 

m c (t e ) = m c (0) - &irm H N o c PDR (r c (0)t e + c c t 2 e /2). (CI) 

Also, 

m c (t e ) = —m H n c {t e )r c 3 {t e ) 3 . (C2) 

Because there is significant expansion, n c (t e ) <C n c0 , we assume that thermal pressure dom- 
inates in both the clump and the PDR flow. From equations (B5) and (C2) 

m c (t e ) = —m H N — r c (t e ) 2 C3 

3 1 + a 
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Equations (CI) and (C3) can be solved for t e , and 

^^-^)\[ , v+ ';z^ + S - 1] (C4) 

At times t > t e , clump expansion is halted, and the clump is now confined by the pressure 
at the base of the PDR flow. The clump slowly loses mass, and shrinks to eventually become 
completely photoevaporated. The time evolution of clump mass and size can be obtained by 
differentiating equation (C3) which also holds for t > t e , with respect to time and equating 
this with the mass loss rate (equation Bll ). Thus, 

3 1 j 

r c (t > t e ) = r c (t e ) - -cpDR—r- — (t - t e ), (C5) 
Z 2z/ z + a 

2(2u 2 + a) fr c (t)\ 2 
m c {t > t e = m c0 1 H — • C6 

The photoevaporation timescale is obtained by setting the clump radius to zero at tpE to 
give 

( 2v 2 + a \ r c (0) + c c t e 

tpE -{-TT^) sc PDR +te ' (C7) 

and t e is given by equation (C4). 



D. Numerical hydrodynamics code 

The fluid equations of motion for the system are solved using a 1-D spherical Lagrangian 
hydrodynamics code. The equations are solved using a finite difference method, with a 
numerical viscosity term, added as a pseudopressure, for accurate handling of shocks and 
discontinuities (Richtmyer & Morton 1967, Bowers & Wilson 1991). An isothermal equation 
of state is used throughout. The accuracy and stability of the scheme was checked with 
standard test problems with known solutions. The computational grid has 1000 equally 
spaced radial zones, except for the central 20%, which was spaced logarithmically in radius. 
This was done to increase spatial resolution at the centre, and it provides a more accurate 
calculation of any shock-induced collapse. 

The FUV field is gradually turned on over a timescale, tpuv, which is varied to accom- 
modate both the impulsive and slow heating cases described in §4. The temperature of an 
outer shell of material is thus raised continuously from T c to Tpdr, reaching its maximum 
value of the PDR temperature at a time tFuv- The temperature T, exponentially drops off 
with the square of the column, N, into the cloud, T{N) =T C + {T PDR - T^e"^'^' 2 . This 
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profile was chosen to closely mimic the T dependence with column predicted by PDR models 
for the values of G and n under consideration (Kaufman et al. 1999). 
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